on Remote Sensing. 
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From the 1970's to the present, scientific investigators of the Office of Naval Research have been using 
environmental satellites such as the geosynchronous satellite shown here. The spin-scan cameras on these 


satellites produce imagery of the earth, its oceans, and clouds. 

ONR has been using remote sensing techniques for many years. Back in the 1950's, ONR studied 
shallow water currents from photographs taken by cameras suspended from helium-filled balloons. Today 
ONR is involved with remote sensing research conducted from manned space crafts equipped with in- 
struments capable of recording oceanographic data. 
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Remote Sensing Research 


by Paul F. Twitchell 
Office of Naval Research, Boston 


ye U.S. Navy has a long history in remote sens- 
ing research and its application to communica- 
tion, navigation, surveillance, weapons detection and 
delivery systems. The Geography Division of the Of- 
fice of Naval Research (ONR) initiated in the 1950’s 
a research program focusing on remotely sensing the 
environment. ONR was a pioneer in this field and 
continues today as an important force in many new 
areas of remote sensing research. In 1962 ONR spon- 
sored the first International Symposium on the 
Remote Sensing of the Environment. In recent years 
these meetings have grown into major conferences 
each resulting in proceedings of several volumes. In 
1973 ONR’s Remote Sensing Program was reviewed 
in Naval Research Reviews’. In this issue are four 
papers focusing on a few facets of remote sensing of 
the ocean. 

For more than forty years, Navy scientists have 
been integrating advances in microwave technology 
into the fleet including shipboard and aircraft fire 
control radars, weapons guidance, and ocean surveil- 
lance from space. In the past twenty years the use of 
radio waves to measure ocean surface dynamics has 
contributed greatly to the understanding of wave 
dynamics. Navy scientists, particularly the late John 
W. Wright of the Naval Research Laboratory, have 
played a significant role in studying phenomena at 
the air-sea interface relating to short-wave modula- 
tion by long waves, the spectral development of 
ocean surface waves, and wave growth predictions. 
Many of the basic scientific questions, however, have 
not been answered; for example, phenomena such as 
nonlinear wave-wave energy transfer and strong non- 
linear wave-current interaction require improved 
theoretical understanding to fully explain radar 
amplitude modulation and to extract wave spectra in- 
formation from radar returns. Professor Erik Mollo- 
Christensen’s article in this issue addresses some of 
these unanswered questions. 

Microwave remote sensing of the ocean uses 
radiometric and altimetry devices. For about a 
decade scanning microwave radiometers have flown 
on U.S. satellites. The Navy has and continues to 


develop techniques for using the ‘‘brightness temper- 
ature’’ as measured by the radiometers. One practical 
application has been the identification and classifica- 
tion of sea ice. Among the information gained from 
radar altimeters on satellites are the locations of large 
ocean eddies and the slope of the ocean surface 
across a major current system. The location and 
speed of currents can be estimated from space using 
ocean slope measurements. 

The infrared sensors that are common today on 
aircraft and satellites are traceable to basic research 
funded in part by the physics and chemistry branches 
of the Office of Naval Research in the 1950’s. Hard- 
ware developments and their applications were ac- 
celerated through the Infrared In’ormation Sym- 
posiums organized and run by OWNR. Special credit 
should be given to two late ONR scientists, T. B. 
Dowd of the Boston Branch Office and L. B. White 
of the Chicago Branch Office for directing these sym- 
posia which provided a forum for the development of 
modern electrooptical weapons systems. 

H. L. Clark and other scientists at the Naval Re- 
search Laboratory developed an early (1942) infrared 
system for remotely sensing the sea surface from 
lighter-than-air ships’. By the 1960’s airborne in- 
frared systems were in common use. Two ONR 
funded investigators, P. M. Saunders at Woods Hole 
Oceanographic Institution and E. D. McAllister at 
Scripps Institute of Oceanography, used airborne in- 
frared systems in their basic investigations of oceanic 
processes. Contemporary ONR-funded efforts that 
employ infrared sensors are reported in this issue by 
Otis B. Brown and Robert Evans of the University of 
Miami and Oscar Huh and Lawrence Rouse of Loui- 
siana State University. 

By using environmental satellites in the early 
1960’s, the Navy pioneered in research which ob- 
tained on a synoptic (1000’s km) and meso (100’s km) 
scale environmental characteristics of the ocean and 
atmosphere. Over the past two decades advances in 
quantifying satellite data have resulted in cost saving 
for data acquisition, and a better understanding of 
the physical processes of the ocean. The Office of 
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Naval Research has sponsored many scientific meet- 
ings and workshops to advance the application of 
remote sensing tools for the study of oceanic proc- 
esses. At the meeting on Physical Oceanography and 
Satellites held at Woodshole in June 1978, Professor 
Verner Suomi concluded with the following 
statement’: 


Remote sensing and data collection from 
satellite platforms show considerable promise 
for improving our understanding of a wide 
range of scales of ocean flows and their in- 
teraction with the atmosphere and underlying 
surface. The data volume from these new 
tools will be immense compared to traditional 
methods of gathering oceanographic data. If 
the full benefits of these new tools and data 
sets are to be realized, new methods are re- 
quired to convert these composite data sets 
into information. The problem is similar to 
that experienced by the meteorological 
satellite operations. We therefore recommend 
that rapid analysis tools using minicomputers 
and interactive video displays be made avail- 
able to oceanographic researchers: further 
that a modest number of these tools be made 
available to individual researchers or to small 
groups rather than through a simple or small 
number of large institutional data processing 
facilities. Our present ability to fully utilize 
these vast data sets is limited. Nevertheless, we 
recommend that as many, and as complete, 
data sets as possible be archived, taking into 
account reasonable cost but at the same time 
being mindful of the future potential value of 
the data for studies of the time-history of large 
scale flows. We stress that any archiving 
scheme will have value only in direct propor- 
tion to provisions for economical and ready 
access to the data of the research and opera- 
tional community. 


The Navy remote sensing program is addressing 
the needs suggested by Professor Suomi and Navy 
environmental support requirements outlined in the 
following table. 

As indicated by Table 1 the Navy environmental 
support requirements span a broad range of time-and 
space-scales demanded by each of the several Navy 
mission areas. Applications of satellite oceano- 
graphic data include support of fleet mission oper- 
ations and support of various research and develop- 
ment programs. 

There exists a Navy Operational Requirement 
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Table 1 
Navy Environmental Support Requirements 


(ADAPTED FROM Noble and Felt) 


MISSION SUPPORT: 
STRIKE, ANTISUBMARINE, 
ANTISURFACE, AMPHIBIOUS, and 
MINE WARFARE 


TIME-SCALES 
HOURS 
DAYS 


(AIR OPERATIONS) 

(TYPHOON FORECASTING, 
ICE OPERATIONS) 

(MISSION PLANNING, SHIP 
ROUTING) 


WEEKS 


SPACE-SCALES 
10m 
1-50km 


(COASTAL PROCESSES) 

(OCEANIC FRONTS, EDDIES, 
CURRENTS) 

(SEVERE STORMS, GLOBAL 
“WEATHER”) 


50-1000km 


APPLICATIONS 
TACTICAL AND STRATEGIC MISSION 
OPERATIONS 
RESEARCH AND DEVELOPMENT 
PROGRAMS 


for Satellite Measurement of Oceanographic Param- 
eters (SMOP OR-WO527-OS). The specifications for 
horizontal resolution, precision, and accuracy of 
measurement of identified parameters are shown in 
Table 2. The dual figures in the statement of reso- 
lution, precision, and accuracy represent estimates of 
minimal acceptable thresholds for operational utility 
when the document was prepared‘. 


To achieve the operational requirements out- 
lined above there is a need for basic and applied 
research on sensors, image enhancement and display, 
data interpretation, integration of satellite data into 
numerical models and other topics. The paper in this 
issue by Thomas Brubecker of Colorado State Uni- 
versity is an example of the efforts to translate in 
near real time the immense volumes of data stream- 
ing down from satellites into a useful presentation 
for an operational person such as an aircraft pilot. In 
a later issue of Naval Research Reviews there will be 
additional articles on improving satellite data trans- 
mission and processing. 

Another area of research is the application of 


Table 2 


Operational Requirements Satellite Measurements of Oceanographic Parameters 


(SMOP OP-W0O527-OS) 


PARAMETER 


HORIZONTAL 
RESOLUTION * 


PRECISION ACCURACY* 


SURFACE WIND 
SPEED 
DIRECTION 


SEA SURFACE 
TEMPERATURE 


SEA VERTICAL TEMP. 
STRUCTURE** 


WAVES 
SIGNIFICANT HEIGHT 
AMPLITUDE COMPONENTS 
WAVELENGTH COMPONENTS 
DIRECTION 


ICE 
COVER 
THICKNESS 
AGE 


ICEBERGS 


10/ /25 km 
10/ /50 km 


10/ /25 km 


10/ /100 km 


10/ /25 km 
1/ /25 km 
1/ /25 km 
1/ /50 km 


0.5/ /25 km 
2/ (50 km 
10/ /50 km 


0.015/ /0.1 km 


2m/sec/ /4m/sec 
+ 10°/ /22.5° 


5%/ 120% 
+5°/ /10° 


0.25/ /0.8°C 0.5/ /1.0°C 


0.25/ /1.0°C 1.0/ /2.0°C 


0.3 m/ /10% 

0.3m/ /0.7m 

5%/ 115% 
10/ /45° 


0.3m/ /10% 

0.3m/ /0.7m 
5%/ 115% 
10°/ /30° 


12%/ 130% 
0.5/ /2m 
6/ /12 mos. 


10%/ 
0.25/ /2m 
6/ /12 mos. 


N/A POSITION 0.5/ /2 km 


* OPERATIONAL TARGET // MINIMUM ACCEPTABLE THRESHOLD 
** DEPTH 300/ /50m; VERTICAL RESOLUTION + 2/ /10m 


satellite imagery of ocean surface characteristics and 
atmospheric profiling into numerical models of 
oceanic and marine boundary layer dynamics. At the 
present time there is a need to find resources to 
continue and expand efforts to optimize the ingestion 
of a satellite data into numerical models used to 
describe and predict operationally important envi- 
ronmental information. 

The remote sensing articles in this issue represent 
only a small fraction of the ongoing remote sensing 
programs of the Navy. In-house remote sensing ac- 
tivities of the Navy include research at Naval Envi- 
ronmental Prediction Research Facility (NEPRF) 
where significant contributions have been made on 
tactical applications of satellite data and the develop- 
ment of a Satellite Processing and Display System 
(SPADS). At the Naval Research Laboratory 
pioneering work on sensor technology is interna- 
tionally recognized. The Naval Surface Weapon’s 


Center, Dahlgren (NSWC/DL) and the Navy Ocean- 
ographic Office (NAVOCEANO) are pursuing 
research using satellite data. At the Naval Ocean 
Research and Development Activity (NORDA) scien- 
tists are developing techniques for effective opera- 
tional use of satellite imagery and investizating 
methods for applying satellite data to numerical 
models. Other in-house efforts related to systeras 
developments include the multispectral satellite uata 
analysis at Naval Ocean Systems Center (NOSC) and 
the application of satellite sea surface temperature 
data at Naval Underwater Systems Center (NUSC). 
Scientists of the West Coast ONR Office have a long 
tradition in fostering new remote sensing cooperative 
efforts with other agencies and various elements of 
the operational Navy. 

Navy scientists and their university colleagues 
funded by the Office of Naval Research and Naval 
Air Systems Command (Air 370) Atmospheric/ 
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Oceanographic Program communicate with each 
other through professional societies. On occasion 
workshops or meetings are organized by ONR to en- 
courage advances in areas of research of interest to 
the Navy. The results of ONR-sponsored activities 
are well promulgated through professional society 
publications or special volumes; recent examples in- 
clude a workshop report in the Bulletin of the Ameri- 
can Meteorological Society’ and a book*®, Navy plan- 
ners and managers draw upon such resources 
through groups such as the Navy Environmental 
Remote Sensing Coordinating Advisory Committee 
(NERSCAC) to introduce research products. For ex- 
ample, a specific research finding could be integrated 
into the Tactical Environmental Support System 
(TESS) development. The Navy operational com- 
munity is presently formulating plans for effective 
and affordable implementation of satellite products 
into global and regional environmental prediction 
systems. 

The author wishes to acknowledge contributions 
from colleagues within the Office of Naval Research, 
the Naval Research Laboratory, and from other 
organizations mentioned throughout the text. Hf 
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apid surveillance of inshore water for bath- 
_.. and underwater obstacles is important 
for a variety of applications. Extensive research 
has been performed in this area during the past 
several years using a pulsed blue-green laser to 
shallow water depth from aircraft. 
Limitations to this system such as turbidity have 
necessitated research on a more versatile or less 
restricting system. 

A new laser technique is being studied for in- 
shore bathymetry and object detection. The 
system is a hybrid system consisting of an infrared 
laser (CO?) transmitter and a highly sensitive 
microphone receiver that is located in the air. 


Continued on page 61 


New Laser Technique for Inshore Bathymetry 


AMPLITUDE - (mV) 


LASER/ACOUSTIC ~=RETURNS 
KEY WEST-AUG 1981 


WATER DEPTH 23ft. 


TIME - (me) 
Figure 1. 
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ACRONYMS 
These Acronyms Appear Frequently in This Issue 

AVHRR Advanced Very High Resolution Radiometer 

CNOC Commander Naval Oceanography Command 

czcs Coastal Zone Color Scanner 

CMSP Defense Meteorological Satellite Program 

ERTS Earth Resources Technology Satellite 

FOV Field of View 

GOSSTCOMP Global Operational Sea Surface Temperature Computation \ 

HIRS High-Resolution Infrared Sounder 

IFOV Instantaneous Field of View 

INDEX Indian Ocean Experiment 

IR Infrared 

IRIS Infrared Interferometer Spectrometer (also Infrared Information 
Symposium) 

LANDSAT Land Satellite (Earth Resources Technology Satellite) 

McIDAS Man-computer Interactive Data Access System 

NASA National Aeronautics and Space Administration 

NAVAIR Naval Air Systems 

NAVOCEANO Navy Oceanographic Office 

NEPRF Naval Environmental Prediction Research Facility 

NERSCAC Navy Environmental Remote Sensing Coordinating Advisory 
Committee 

NESS National Earth Satellite Service 

NIMBUS Polarorbiting Meteorological Satellite 

NOAA National Oceanic and Atmospheric Administration 

NOAA-5, NOAA-6, NOAA-7, etc. Operational Polarorbiting Meteorological Satellites 

NOMSS Naval Oceanographic and Meteorological Support System 

NORDA Naval Oceanographic Research and Development Activity 

NOSC Naval Ocean Systems Center 

NRL Naval Research Laboratory 

NSWC Naval Surface Weapons Center 

NVSC Naval Underwater Systems Center 

RSMAS Rosentiel School of Marine and Atmospheric Science 

SAR Synthetic Aperature Radar 

SASS SEASAT-A Scatterometer System 

SCAMS Scanning Microwave Radiometer 

SEASAT Sea Satellite (a satellite designed to measure oceanic . 
characteristics, it operated for a short time in 1978) 

SEBEX Surface Expression of Bathemetry Experiment 

SLAR Side Looking Airborne Radar 

SMMR Scanning Multifrequency Microwave Radiometer 

SST Sea Surface Temperature 

TESS Tactical Environmental Support System 

THIR Temperature, Humidity, Infrared Radiometer 

VAS Vertical Atmospheric Sounder 

VTPR Vertical Temperature Profile Radiometer 
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greater importance to the understanding of oceanic 
variability than the quasi-quantitative values of the 
data. 

Satellite-based remote sensing of the environ- 
ment first reached public notice with display of color 
earth pictures taken from the Mercury manned 
spacecraft. In fact scientific sounding rockets and 
military passive sensors were routinely mapping the 
oceans before this time'?*. Multi-spectral color pic- 
tures captured public and scientific fancy. However, 
it was to be some 15 years later with the launch of the 
NIMBUS-7 Coastal Zone Color Scanner before these 
promises could start to be realized over oceanic 
areas; the Earth Resources Technology Satellite, later 
called LANDSAT, provided quasi-operational land 
remote sensing in the early 1970s. 

Application of infrared remote sensing to the 
study of meteorological variability was quickly ap- 
preciated by analysts and an operational military 
program was initiated in 1963. This program was 
managed by the U.S. Air Force and concentrated on 
determination of world-wide cloudiness, cloud-top 
temperatures and tactical applications for aircraft 
logistical use. It was not focussed on the study of 
oceanic surface thermal variability. The past decade 
has seen active Navy participation in these efforts, 
now called Defense Meteorological Satellite Pro- 
gram, with a gradual orientation of the system 
toward improving ocean thermal-imaging capabil- 
ities. 

In this report we will discuss the visible/infrared 
sensing problem, some applications and conclusions. 
The sensing problem logically splits into two sec- 
tions: technology to acquire the data, ie., sensors, 
and algorithms or processing paths for the data once 
acquired. We will not discuss the sensor problem ex- 
cept as needed for algorithm development. 


Sensor Processing 


Sensor data processing can be divided into the 
following 4 categories: 

1) ingest, calibration, and remapping 

2) navigation 

3) calibrated image retrieval, magnification, 

geometrical correction, and hard copy output 

4) interactive analysis with ground truth data, 

contour plots, etc. 

These four classes of activity are necessary for 
any satellite data analysis procedure; similar systems 
have been developed for LANDSAT, SMS, DMSP 
and other geophysical sensing satellite programs. In 
fact, as will be noted, parts of the procedures to be 
described are adapted from other efforts. The pres- 


ent system is oriented towards polar-orbiting, rotat- 
ing scanner, small acceptance-angle sensors. Figure 1 
summarizes the overall procedure. 


Data Ingest and Calibration 


A qualitative appreciation for passive remote 
sensing measurement techniques is necessary to 
understand properly the data ingest procedure 
design. Consequently we will give a brief introduc- 
tion to visible/infrared remote sensing technique 
prior to discussion of data ingest procedures. 

Remote sensing systems make radiometric mea- 
surements in specified electromagnetic frequency 
bands of energy emitted or reflected by the earth. 
The absorption spectrum for radiation passing 
through the atmosphere to space from the sea surface 
is shown in Figure 2. There are three areas of this 
spectrum (‘‘windows’’) in which we have interest: the 
visible band between 0.3 and 0.7 micrometers, and 
infrared bands between 3.7 and 4.2 micrometers and 
10.0 and 12.0 micrometers. These windows all have 
the characteristic of having relatively low absorption 
for radiation; thus they yield relatively clear views of 


TYPICAL SATELLITE DATA PROCESSING PROCEDURE 


FoR VHRR, SR, DMSP, GOES 


ACQUISITION 


REFORMAT FROM SOURCE MEDIA TO INTERNAL FILE 


CALIBRATION 


1. PRELIMINARY VALIDATION 2. SENSOR CHARACTERISTICS 


3, ELECTRONIC CHARACTERISTICS 


NAVIGATION/RECTIFICATION 


1. EQUAL AREA PROJECTION 2. COORDINATE CONVERSION 


3. LANDMARK 4, DESKEW IMAGE 


DATA SASE COMPILATION 


OEVELOP IMAGE FILES FOR DESIRED COMPOSITION, RESOLUTION 


ENHANCEMENT FILTER COMPOSITE PLOT/ UNIVAC TRANSFER 


PRINT 


Schematic view of the satellite data processing pro- 
cedure. 


Figure 1. 
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the earth’s surface. Scanning radiometers operating 
in these windows produce two-dimensional images of 
the sensed outgoing energy. A point value at the 
spacecraft for a scan-line element pair (pixel) is 
sensed as either a voltage or current value propor- 
tional to the incident energy. This analog value is 
digitized and telemetered to earth, or telemetered and 
digitized on earth. TIROS-N is an example of the 
former*, and NOAA-S5 of the latter’. Final digital 
count values are transcribed into computer-compati- 
ble format on earth. Possible error sources in this 
system occur during the digitization and transmission 
steps. NASA, NOAA and DMSP spacecraft-data 
streams contain known references by which telemetry 
and digitizer errors can be identified and corrected. 
Sensor calibration is subject to drift as well, due to 
changes in environmental factors aboard the space- 
craft. This three-step process can be described func- 
tionally: 


L = L(C(V(A))) (1) 


where L is a radiance, an energy measure, C a drift- 
voltage transfer function, V a voltage transfer func- 
tion, and A the telemetered signal leVel. In the case of 
infrared radiometers the final temperature is a non- 
linear function of the radiance L, it is not directly 
sensed. Seawater can be characterized as a grey body 
emitter in these two infrared windows with an emis- 
sivity, € of .98. To convert the energy value to 
temperature, a blackbody transform is applied to L, 
i.e. 


T=T{(L/e) 2. 
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Figure 2. Atmospheric transparency for sea surface light as seen 
from a satellite. Computed for a U.S. standard at- 
mosphere utilizing the LOWTRAN 4 code (from Selby 
et al., 1978). 


where T is an equivalent blackbody temperature 
(Planck) function. These transforms have been 
documented by Koffler®. 

T is a radiometrically correct temperature at the 
top of the atmosphere. Two principal effects must be 
taken into account: variation in atmospheric absorp- 
tion as a function of scanner angle (path length) and 
humidity. Path length absorption is strongly con- 
trolled by the atmospheric water vapor since the ver- 
tical distribution of water vapor is not constant. 
Scientists have examined these effects and shown that 
the sensed temperature can be up to 10C in error for 
very humid maritime atmospheres at large scan 
angles.’’* Two methods are presently available in our 
laboratory for calibration. The first is a zero-order 
method based on climatological mean humidity pro- 
files. This was developed by National Earth Satellite 
Service and is described by Smith’. A tempera- 
ture correction matrix is derived by regression as a 
function of scan angle and sensed surface tempera- 
ture. Sensed surface temperature and total water 
vapor concentration are found to be highly cor- 
related. Scenes with homogeneous moisture fields are 
readily analyzed by this technique which produces 
semi-quantitative two-dimensional surface tempera- 
ture maps. Functionally, 


T = H(a, T) (3) 


with a the scanner angle and T calculated from equa- 
tion 2. This correction, though based on climatic 
mean moisture profiles, is ideal for quick-look and 
image sorting applications. 
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The second method which provides a better cor- 
rection for atmospheric absorption is based on multi- 
spectral sensing in the water vapor bands. The new 
AVHRR radiometers have two or three channels in 
the infrared: 3.7y-4.2u, 10y-11p and 11-12, which 
are channels 3, 4 and 5 respectively. The integrated 
atmospheric transmissivity over each of these three 
channels differs, and consequently algorithms can be 
constructed which depend on the difference of the 
perceived surface temperature in these channels. The 
simplest algorithm assumes that the atmosphere is 
optically thin enough that the perceived surface 
temperature difference between any of the channels 
can be parameterized as a simple function of the bulk 
atmospheric water vapor seen by the satellite sensor. 
Functionally a formula of the following form is ap- 
plied: 


+ (Fs = T,) (3a) 


where the T,’s are temperature retrievals in various 
channels for a given location and the a;’s give the 
parameterized correction. The a;’s can be derived 
from first principles theoretically'® or can be derived 
empirically from good composite sets of surface and 
satellite observations.'' Both of these methods give 
better control on the temperature accuracy at a 
sacrifice in processing throughput. 

Visible measurements suffer from similar prob- 
lems although they arise from different sources. An 
analyst rarely desires a radiance measurement at the 
top of the atmosphere, but rather wishes to know the 
upwelled radiance at the sea surface or the mean 
chlorophyll concentration or sediment load in sur- 
face waters. Whereas emitted thermal radiation is an 
intrinsic property of the medium, upwelled visible 
light arises from back-scattered solar radiation. Con- 
sequently, one must calculate the amounts of solar 
radiation reaching the surface directly and scattered 
from the atmosphere, the resultant that is trans- 
mitted into the sea, that which is back-scattered by 
the sea and suspended materials and is transmitted 
through the surface and back through the atmo- 
sphere to the satellite. The atmosphere, ocean and 
suspended material selectively affect the down and 
upwelled radiances as they propagate from the 
sources to the satellite sensor. It is a problem which 
has recently been straightforwardly parameterized. 
Atmospheric correction algorithms have been pro- 
posed.'?:'*:'4'5 The method of correcting for these 
effects'? and with the clear water modifications'® is 
summarized in the following. 

Radiance which reaches the sensor L,(A) is com- 
posed of contributions from atmospheric Rayleigh 


scattering L,(A), aerosol scattering L,(A) and finally 
the desired contribution L,(A) from the water. These 
terms combine as: 


where t(A) is the diffuse transmittance of the atmo- 
sphere. This formulation assumes that the sensor is 
not viewing specularly reflected sunlight (glitter) and 
that the fresnel reflection from the sea surface is in- 
cluded in the Rayleigh/aerosol terms. The Rayleigh- 
scattered light is directly calculable from theory. 
Hence, if L(A) is known, then L,(A) is derivable. The 
initial formulation of the correction process uses this 
approach.'’ If it is known that the aerosol size- 
refractive index distribution is horizontally homoge- 
neous over the region of interest, then this algorithm 
quickly generalizes. 
This implies that the aerosol selectivity: 


= Ly 2) / L, Os) (5) 


is constant over the region of interest. 

In offshore areas, such as the Sargasso Sea or 
the Gulf Stream, it is well known that there is very lit- 
tle terrageneous material in surface waters; Figure 3 
shows that upwelled visible radiance in the red (670 
nanometer channel for NIMBUS-7) is quite small as 
a consequence. One thus can assert that whatever is 
seen in this part of the spectrum must arise from the 
atmosphere. That is, we assume 


tA) Ly) = 0 (6) 


Then in a straightforward way one finds that 
LyQ) = {L,Q) - 
SA, Ap) — Lp) (7) 


Alternatively, S(A, \,) can be computed with high ac- 
curacy if L,(A) is known for some area in the scene. 
Gordon and Clark'® show that L,(A) for Case 1 
waters can be calculated with just a knowledge of the 
solar zenith angle and that operationally these waters 
are easily identified.'* Thus, in this method, the pro- 
cedure is to find a ‘‘clear water’’ area, compute L,,(\), 
derive S(A, \)), and proceed to correct the jmage. A 
case study is presented later in this report. 


Geometrical Correction 


Equally as important as radiometric correction 
to the final data product is conversion of scanner 
elements to some easily manipulated aspect ratio 
(metric). Some satellite systems, e.g. DMSP, ac- 
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Figure 3. Upwelled light spectral sensitivity to suspended mate- 
rials (from Gordon and Clark, 1980).*’ 


complish a portion of this in hardware by control of 
scanner or display write-rate'’®. NOAA and most 
other polar orbiting satellites have fixed-rate scan- 
ners and are operated with orbits that have inclina- 
tions differing from 90 degrees*. Thus, these in- 
struments geometrically distort the pixel aspect ratio 
at larger scanner angles. Other distortions due to 
satellite vehicle pitch, roll and yaw may also be pre- 
sent. Pitch, roll and yaw have much more compli- 
cated interactions with the scan-line aspect ratio. 
Variations in these parameters are deduced by study 
of the differences in high precision navigation 
generated and observed landmark positions. 
Legeckis and Pritchard*® describe an algorithm 
for geometrical correction of scan line data. Al- 
though similar methodology is used, latitude, time 
and height-dependent functionals are not based on 
mean estimates, but carry the full dependence of each 
parameter. Estimates of these parameters are avail- 
able from an ephemeris-driven, high precision navi- 
gation algorithm that will be described in a latter sec- 
tion. The geometrical correction produces scan lines 
containing pixels with a prescribed aspect ratio. 
Functionally a unidirectional mapping from scan- 
element to new constant aspect ratio scan-element 
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pixels is defined: 


F = F(s, e, , h) (8) 

where s is scan line, e element, @ latitude and h 
satellite height. Figure 4 sketches the algorithm. To 
first order, there are two effects to be corrected: pixel 
distortion along a scan line and relative scan line mo- 
tion. The pixel distortion is taken to be discrete, i.e., 
old samples are mapped to the nearest whole new 
sample. As shown in Figure 5, the vector contains 
varying repetition factors ‘‘n’’, which increase as a 
function of distance from the nadir point. A given 
old sample is mapped into ‘‘n’’ new samples. A single 
run-length encoded transform vector is constructed 
based at the nadir point, since, to first order, the cor- 
rection is symmetric about the nadir track. This 
transform is quite fast due to run-length encoding of 
the mapping. The second distorting factor, earth 
rotation, to first order, manifests itself as a line-shift 
correction. These are either line slips or advances 
relative to the previous line dependent on satellite 
nadir track sense, scanner mode and hemisphere. 

Data sets that result from transforms 1-8 com- 
prise remapped scan-line sets which are stored on 
disk and/or tape. These fixed-format, calibrated, 


11 


| 
: 
) : 10 
1 
4 we, 
x 
~ 
5 
| 
+ 


IMAGE RECTIFICATION 
EQUAL AREA MAPPING 


37 3/2/2/1/2) 2/3 


EARTH ROTATION CORRECTION 


- 


Figure 4. Pictorial representation of scan line and earth rotation 


remapping corrections. 


constant-aspect-ratio data sets constructed from the 
polar orbiter raw data, are geometrically and radio- 
metrically similar to geosynchronous spin-scan and 
other small acceptance angle scanners. All subse- 
quent analyses utilize this data-base store for input 
which simplifies the succeeding levels of analysis 
algorithms. 


High Precision Navigation 


Orbiting vehicles can easily be engineered with 
accurate time bases, but it is more difficult to main- 
tain precise vehicle position and attitude informa- 
tion. Older sensing systems, of which NOAA-S is a 
prime example, have marginal attitude sensors (+ .5 
degree stabilization)*. Newer vehicles such as TIROS- 
N‘’, DMSP-BLOCK 5D,'* have accurate attitude 
(+ 1° stabilization) and temporal (+ 50 msec.) deter- 
mination systems. The TIROS-N time-keeping ac- 
curacy is primarily determined by ground-station 
operational constraints on spacecraft clock update 


12 


frequency. Since a one-second error translates into 
an along-track position error of 6.6 km., this is a 
tight constraint on update accuracy. We observe 
(1979-1980) variations as large as two seconds (13 km 
along track error) in clock accuracy. 

In all cases there has been little done to provide 
useable onboard navigation-linked sensor process- 
ing. An ephemeris-driven algorithm,' has been con- 
structed to provide latitude, longitude and height in- 
formation for low-precision mapping. A Keplerian 
model is utilized with second and fourth order grav- 
ity terms; second order secular perturbation correc- 
tions on the period, mean anomaly and ascending 
node are applied in addition. Functionally the posi- 
tion: 

Py 6h = (Ephemeris) (9) 
is derived. The ephemeris consists of the orbital 
elements: longitude of ascending node, angle of or- 
bital inclination, argument of perigee, semi-major 
axis, eccentricity, and universal time. \’ is the 
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Figure 5. Cruise track and derived temperature field from R/V 
Columbus Iselin Cruise CI7704 for April 10-12, 1977 (Courtesy Dr. T. Lee). 
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longitude, ¢ latitude, h height. This procedure is 
well known and is lucidly documented”’. Tests of this 
high-precision navigation algorithm, P ,.,,, demon- 
strate that satellite position can be predicted on the 
basis of daily mean orbital elements (ephemerides) to 
+ 1 km. 


Precision Navigation/Mapping 


The high precision of the ephemeris-driven 
navigation cannot be fully utilized to establish precise 
scan, element location due to spacecraft clock drift 
(c-13 km induced error) and attitude control(+ 4 km 
induced error in roll, pitch, and yaw). Quantification 
of these errors requires generation of a precision 
polar-orbit navigation model. Philips and Smith”? 
describe a landmarking (control point) approach for 
navigation of geosynchronous satellites. Knowledge 
of the scanner geometry permits computation of cor- 
rections to the ephemerides by comparison of known 
earth location positions with model-predicted posi- 
tions. Daily update of the model then allows ‘‘predic- 
tive’’ navigation of the spin-scan radiometer data. 
This approach is not directly applicable to the low 
altitude, several-pass-per-day polar-orbiter situation 
due to the presence of non-secular perturbations. 

A navigation method for polar orbiters has been 
developed in our laboratory.** In this method, the 
nadir track and the scan lines are considered to define 
great circles on a spherical earth such that the any 
scan line great circle is always perpendicular to the 
nadir track great circle. (This is equivalent to the ap- 
proach of Ruff and Gruber,*’ however, the trans- 
form analysis is performed quite differently.) The 
procedure to establish navigation parameters for a 
satellite pass consists of the following steps: First, 
from the universal time of the initial scan line, the 
high precision navigation model is used to calculate 
the earth location of the nadir pixel. Second, land- 
marks (or continental outlines) are selected with an 
emphasis on maximizing their spatial spread with 
respect to the nadir track. This is done to minimize 
the entry of systematic error in the transform equa- 
tions. Third, corrections from the outlines can be 
computed for the coordinates of the initial scan line, 
the time correction and mean pitch, roll and yaw of 
the satellite. 

Analytical transforms are generated by this tech- 
nique which relate constant-aspect-ratio pixel scan- 
line elements to latitude, longitude and vice-versa. 
These transforms form a basis for derivation of 
mapped projections of scanner data. 

Performance of detailed analyses of satellite 
derived ocean skin temperature also requires a 


humanized interface for data extraction, enhance- 
ment and display. An intelligent extraction algorithm 
founded on the precision navigation transforms 
yields an easily useable method to display an image 
centered at some earth/scanner coordinates with a 
given resolution. This resultant data set (or the 
original) with an appropriate transform function can 
then be projected into various coordinate systems 
(transverse mercator, polar stereographic, etc.) 
which are stored on disks and are available for recall 
at the analyst’s convenience. 


Case Studies 


In the first part of this report we summarized the 
radiative transfer, navigation and data bases needed 
for satellite remote sensing of visible and infrared ra- 
diance from the ocean. Now we will describe several 
case studies which illustrate that application of the 
above methodology yields quantitative and qualita- 
tive improvement of our understanding of the ocean. 


NOAA-5 Gulf Stream Location Observation 
in the South Atlantic Bight 


During early April, 1977, Dr. T. Lee (University 
of Miami/Rosenstiel School of Marine and At- 
mospheric Science) mapped the thermal structure of 
the Gulf Stream/shelf water interface in the Georgia 
Bight.’ A surface thermal field derived from these 
measurements is shown in Figure 5. Simultaneously 
the ocean skin temperature was observed from space 
by the NOAA-5 Very High Resolution Infrared 
Radio-meter. These data were processed by the 
afore-mentioned procedure and a 3 km gridded 
equal-area product was generated for analysis pur- 
poses. This result is illustrated in Figure 6. All con- 
tinental outlines and political boundaries have been 
inserted from a national data base in terms of earth 
location pairs via the earth-scanner mapping 
transform. Retrieval of the temperature values along 
the ship track is accomplished with the inverse mapp- 
ing transform. The ship track location and cor- 
responding temperatures were digitized as triplets of 
latitude, longitude and temperature from Figure 5 
with a graphics tablet. Figure 7 shows the highlighted 
ship track off the Georgia coast. Magnifying the area 
of interest about the ship track (Figure 8), we notice 
the strong gradients in thermal field at the Gulf 
Stream/shelf water interface. Visual comparison of 
Figures 5 and 7 illustrates the fact that the satellite 
data with its inherently higher spatial definition can 
yield a good surface thermal map in situations with 
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Figure 6. NOAA 5 VHRR equal-area processed thermal map for 


10 April 1977. 


low to moderate atmospheric water vapor (< 1 
gm/cm?). [A mean airmass correction of 1.25 C is 
subtracted from the satellite observations.] The 
temperature along the ship track comprises a one 
dimensional series as a function of cumulative along- 
track distance. Frontal definition and temperature 
agreement between both methods is quite good. 
Quantitatively we find the standard deviations of the 
ship and satellite track-line series to be 2.1C and 
2.2C, respectively, while the rms difference of the 
two series is 1.2C. These values are calculated for a 
set of several thousand observation pairs. The near 
equivalence of the ship and satellite standard devia- 
tions implies there is little second order-moment bias, 
or conversely the satellite shows the same biases as 
the ship. A mean square difference of 1.2C is the ex- 
pected difference at + 1 standard deviation. This is 
comparable to estimates NOAA/NESS has made for 
other temperature products.?’ We believe the main 
contributor to this difference is the inherentiy dif- 
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ferent sampling scales of the two methods. The mean 
square and mean differences are extraordinary in 
that the comparison values are representative of total 
error, that is errors arising from navigation inac- 
curacies, time and space sampling mismatch, water 
vapor contamination, radiometer calibration, etc. 
This verification area was chosen for two reasons: 
high quality ship and satellite data were available 
nearly concurrently and strong thermal features with 
small scales were present which would magnify 
navigation, mapping, and calibration errors. Thus 
we believe that this comparison is indicative of the 
procedure’s utility in relatively clear, homogeneous 
atmospheric water vapor situations. A spatially in- 
homogeneous moisture field would render this 
analysis less valid. 
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NIMBUS-7 CZCS Observations of Color 
Variability in the Middle Atlantic Bight 


A recent work by Gordon et al*® utilizes CZCS 
(Coastal Zone Color Scanner) observations to map 
the surface chlorophyll structure of a warm-core ring 
off the Northeast coast of the United States. Warm- 
core rings occur when extended meanders of the Gulf 
Stream pinch off on the cold-wall (north side) and 
form a closed anticyclonic rotating vortex of Gulf 
Stream and Sargasso Sea waters. This vortex of Gulf 
Stream and Sargasso Sea waters is imbedded in the 
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Figure 7. Magnification of Figure 6 about the ship track. 


surrounding slope or shelf waters and is thought to 
isolate an ecosystem at some point in its life cycle. 
Rings have horizontal dimensions of 100-200km and 
vertical extents equal to or greater than 1000m.’°*° 
We show in this study that mean chlorophyll concen- 
tration in surface waters can be observed to 35% 
utilizing satellite method alone. 

Figure 9 is an overview of the shelf and slope 
water chlorophyll concentration off the northeast 
United States coast for 15 June 1981. The readers at- 
tention is directed to several prominent features: the 
low chlorophyll values (lighter grays) in a warm core 
ring located at 39° N, 69° W, (catalogued as 79-B by 
Fitzgerald and Chamberlin’') the low and high chlo- 
rophyll streamers wrapping into the ring, their 
length, and the overall decrease of chlorophyll as one 
travels north of the Gulf Stream and away from the 
coast. 


OBSERVATION No. 


INTERCOMPARISON OF SHIP VS. SATELLITE SST FOR 04/10/77 


Figure 8. Satellite versus ship derived temperature intercomparison. 
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The fact that this warm core ring has such a 
strong chlorophyll signature implies that it can be 
tracked as either a color or an infrared anomaly. 
Qualitatively the color structure raises several in- 
teresting theoretical questions: How is it that the sur- 
face chlorophyll streamers cross the ring-front and 
spiral in the central regions of the ring? What is the 
mechanism underlying propagation of ring effects 
(through the streamers) to distances of one to two 
ring-diameters away? Are there differences in the 
thermal and color signatures of a warm-core ring? 
Does a regional water mass analysis support the 
large-scale chlorophyll surface pattern we observe? 
Comprehensive delineation of the first two effects 
awaits 1982 warm core ring cruise work, the latter 
question can be addressed qualitatively by looking at 
the slope water mass distributions in the area. 
McLellan’s analysis*? of the water mass areal distri- 
butions in this area agrees with the present observa- 
tions in that he finds shelf-type waters in the Bight 


and up into the Gulf of Maine and cooler, lower ox- 
ygen waters of Labrador origin offshore. This jux- 
taposition of water masses yields the strong horizon- 
tal chlorophyll contrasts seen in the present imagery. 

Figure 10 is an enlarged view of Figure 9 
centered about the ring. Streamer delineation and 
local chlorophyll minimum in the ring interior are 
seen. Figure 11 is an enlarged view of the ring as seen 
several hours earlier in TIROS-N AVHRR infrared 
imagery. Ring frontal outlines are consistent between 
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the different sensing methods, however delineation 
of cross boundary flows is quite a bit more difficult 
in the infrared. The local minimum in central waters 
seen in the chlorophyll is mirrored in the thermal im- 
agery as a temperature maximum. One is left with the 
impression that although there is coherency between 
the two methods, color is a more sensitive indicator 
for study of cross-boundary surface flows. 
Fine-scale features may be atmospheric arti- 
facts. One must know which features are atmo- 
spheric and/or oceanic in origin. A cruise in support 
of the NIMBUS CZCS validation made a section 
through this ring at 0000 EDT, 10 June 1979. NIM- 
BUS orbits 3157 and 3171 bracket this ship coverage 
occurring approximately 12 hours before and after, 
respectively, ship passage through ring center. 
Pigment concentration was measured along the 
ship track by using a Turner Designs flow-through 
fluorometer (Model 10) to monitor the in-vivo fluo- 
rescence of water drawn continuously from about 


Figure 9. NIMBUS-7 CZCS derived chlorophyll image for 15 
June 1979, orbit 3226. Lighter shades denote lower 
values, darker higher (from Gordon et al., 1982). 


1.5m below the surface. Periodically, discrete 


samples were drawn from the same source and ex- 
tractions performed to determine fluorometrically 
the concentrations of chlorophyll and the associated 
phaeopigments, using the technique described by 
Yentsch and Menzel**? with the modifications of 
Holm-Hansen et al**. These discrete determinations 
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Figure 10. 15 June 1979 CZCS-derived chlorophyll image magni- 
fied about a warm core ring found at 39° N., 69° W. 
Line indicates ship track. 


were then used to relate the continuous in-vivo fluo- 
rescence measurements to the pigment concentration. 
CZCS-derived pigment concentration (heavy line) is 
compared with these continuous determinations in 
Figure 12. This track is from south to north, starting 
in the Gulf Stream, passing through the slope water, 
the warm-core ring and finally ending on the shelf. 
These distinct water types are apparent in Figure 12. 
The cold wall of the Gulf Stream is at about 220km 
on the axis, and the ring extends from 300 to 350km. 
Satellite determined pigment concentration is split 
between orbit 3157 from 320 to 440km and orbit 3171 
from 420 to 600km. Gaps in the satellite-derived 
values are due to clouds. Agreement between the two 
methods is apparent. 

Examination of warm-core rings is an ideal ap- 
plication for remote sensing since they have signa- 
tures in the color and infrared. Technically they are 
easy to study as well since the Gulf Stream offers a 
clear-water area for CZCS calibration and slowly 
varying thermal signal for infrared observations. 


Figure 11. TIROS-N thermal infrared image for 15 June 1979 
magnified about the warm core ring shown in the 
CZCS chlorophyll image, Figure 10. Warm temper- 
atures are lighter, cooler darker. Line denotes ship 
track. 


TIROS-N/NOAA-6 Observations of 
Somali Current Response to 
the Summer Monsoon 


Our purpose in this section is to give a summary 
of surface temperature variability near the East 
African Coast as observed by satellite and ship dur- 
ing the mature phase of the summer monsoons of 
1976-1979. Sea surface temperature (SST) pattern 
changes have been shown to be strongly correlated 
with changes in the surface-current field****; as such 
they yield a sequence of observations of the East 
African current-system variability during the summer 
monsoon. 

Monsoonal forcing has two major large-scale ef- 
fects on the northwestern Indian Ocean: Somali Cur- 
rent spinup and maintenance, and cooling of the up- 
per mixed layer. Delineation of these effects on the 
ocean is complicated by the presence of large scale 
eddies near the ocean boundary (which maybe the 
Somali Current) and in the ocean interior®’, intensive 
upwelling near the boundaries**, and other large scale 
disturbances. 

Historically physical oceanographers have con- 
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Figure 12. Comparison of ship and satellite derived surface pig- 
ment concentrations for NIMBUS-7 orbits 3157 and 


3171, and R/V Athena II in-situ observations. Ship 
observations are denoted by the light line, satellite 
observations by the heavier solid and dashed lines, or- 
bits 3157 and 3171 respectively (from Gordon et al., 


1982). 


sidered the Somali Current to be a seasonally revers- 
ing western boundary current. Source waters from 
south of the equator flow north along the East 
African Coast, separating and flowing into the inte- 
rior as part of the ‘‘Great Whirl’’ between 10° and 
12 °N?°4°-4!-42 the Great Whirl being a large anti- 
cyclonic eddy which is observed off the coast of 
Somalia during the summer monsoon. It has dimen- 
sions of 500 km (parallel to coast) by 400 km (perpen- 
dicular to the coast.*? Bruce and Brown find several 
instances (1970, 1976, 1979) when the Somali Current 
is observed to separate from the coast south of the 
Great Whirl between latitudes 2° N-5° N.**:*5 It was 
assumed that there were ‘‘single gyre years’? and 
“*two gyre years.’’ In 1979 Brown found both cases at 
different stages of the summer monsoon.** 

Figure 13 illustrates the ship derived surface 
thermal fields as observed during the mature phase of 
the summer monsoon in 1979**. In the surface- 
current field at this time offshore flow was found at 
2°N, together with an anti-cyclonic gyre between 5° 
N and 10° N with offshore flow from 10° N to 12° 
N.** Surface current and thermal fields are strongly 

‘coupled as might be expected from earlier observa- 
tions**. Separation regions are expressed in the SST 
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as cold wedge shaped areas. 

Figure 14 presents a time series of frontal posi- 
tions as observed during the period June to October, 
1979, by satellite. The southern frontal region starts 
to move poleward, slows and then quickly translates 
to a position ~1° south of the original northerly 
separation location. During the same time interval 
the northerly region (dashed area) moves slowly 
poleward to the region of Socotra and thence more 
quickly southward, coalescing with the southern 
region south of its earlier position near Ras Hafun. 
Ship observations made during the latter phases of 
this process near Ras Hafun find fresher, cooler 
water in the upper ~100 m. overlaying warmer, 
more saline waters in the northwest quadrant of the 
Great Whirl. That is, surface waters are characteris- 
tic of Somali Current origin while deeper waters are 
compatible with Arabian Sea sources (John Bruce, 
private communication, Woods Hole Oceanographic 
Institution). 

Analysis of the 1979 results stimulated interest in 
reprocessing of earlier satellite data archived during 
1976-1978. Figure 15 illustrates observed frontal time 
series for 1976. Southern frontal translation occurs 
later and is much slower than the 1979 case. The 
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Figure 13. Time sequence of ship observed sea surface temper- 
ature fields for oceanic response to the southwest mon- 
soon, 1979. la: 16 April-10 May, 1b: 11 May-22 May, 
Ic: 23 May-12 June, Id: 22 June-31 July, le: 18 
August-3 September. (from Brown et al, 1980). 


multi-phase nature of the transition is evident, al- 
though the absolute changes in translation rate are 
much less. 

Satellite imaging of the area in 1977 was inade- 
quate for resolution of weaker thermal differences 
due to spacecraft malfunctions and related problems. 
Imagery on hand neither supports nor denies the 
presence of the dual separation region in this year. 
Research ship observations of surface temperature 
are available for 1977 from the R/V Shokals’kii 
which was deploying INDEX moorings off Somalia. 
A temperature map for early June is presented in 
Figure 16 (Donald Johnson, private communication, 
Old Dominion University, Norfolk). Temperatures 


are not as low as in the other three years, but a wedge 
is apparent at 6 N. 

Satellite-derived frontal positions during sum- 
mer 1978 are illustrated in Figure 17. Fronts occur at 
the same times as in 1979, however, the northern 
frontal region passes through the Socotra topo- 
graphic arc and eventually dissipates in the western 
Arabian Sea. Southern wedge translation rates are in- 
termediate between those observed in 1976 and 1979. 

Satellite-derived frontal patterns suggest that 
dual separation regions (double gyre configuration) 
occur in three of the four years studied. This inter- 
pretation is at variance with classical views of Somali 
Current circulation. 
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Figure 14. Time sequence of satellite observed frontal locations 
for northern and southern frontal wedges for year 
1979. The northern wedge front is hatched (from 
Evans and Brown, 1981).** 
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Figure 15. Time sequence of satellite observed frontal locations 
for northern and southern frontal wedges for year 
1976. The northern wedge front is hatched (from 
Evans and Brown, 1981). 
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Figure 16. Surface temperature map derived from measurements 
made onboard the R/V Shokals’kii during May-June, 
1977 off Somalia. (from Evans and Brown, 1981, 
courtesy Dr. D. Johnson). 


These observations show coalescence of the 
Somali Current system is abrupt and can be thought 
to consist of four phases. During phase | the 
southern separation region moves north from its 
onset location during the mature phase of the sum- 
mer monsoon to ~6° N. Simultaneously the north- 
ern frontal region translates poleward to the region 
of Socotra. In phase II the system is in quasi-steady 
state for a period of several days to a week. The third 
phase begins with reversal of northern frontal propa- 
gation in 1976 and 1979; in 1978 phase III is initiated 
when the northern frontal zone breaks through the 
topographic arc connecting Socotra and East Africa. 
The southern region moves more quickly poleward 
than in phase I, while the northern zone translates 
equatorward (1976, 1979); in 1978 the northern zone 
translates into the Gulf of Aden and dissipates. 
Finally coalescence occurs during phase IV in 1976 
and 1979. This phase lasts for periods of days to 
weeks. 

From the data presented for the period 1976- 
1979 we assert that this ‘‘anomalous’’ separation is 
confirmed in 3 of the 4 years. More interestingly 
these satellite data give us a first look at the temporal 
evolution of the Somali Current from the anomalous 
separation case to the more historically expected cir- 
culation. It may well be that the single or dual gyre 
nature historically observed is a function of the time 
when observations were made rather than being in- 
dicative of interannual differences in oceanic 
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Figure 17. Time sequence of satellite observed frontal locations 
for northern and southern frontal wedges for year 
1978. The northern wedge front is hatched (from 
Evans and Brown, 1981). 


response. 

Spinup model results suggest that, at least in 
some years, the Somali Current can separate from 
the coast south of the Great Whirl at 5° N. This 
variability may be due to the strength of the mon- 
soon. Cox has shown the generation and stall char- 
acteristics of the eddy circulations are strongly 
dependent on the the wind stress.** If one assumes 
that there is a transition from a stalled-eddy state to a 
poleward-translating state with diminution of the 
wind stress, then the Cox (1979) model may also be 
relevant to the phase I observations. 


Conclusions 


In this report we give an overview of the sensing 
data flow problem, solutions and several applica- 
tions. We have shown that it is possible to do quan- 
titative visible and infrared remote sensing and apply 
it to an improved understanding of the ocean. How- 
ever, an important facet of these applications is that 
it is the patterns of a given parameter which yield the 
most valuable information, not the quantitative 
results. For example, one asks: Where is the Gulf 
Stream? not, How warm is it? or, Where is a warm 
core ring? and, What is its size? Does the ring have 
high or low chlorophyll values? What is the pattern 


of SST response to the summer monsoon? This is not 
to say that quantitative results are unnecessary; they 
directly give us an appreciation for atmospheric con- 
tamination of retrieval and enable us to judge which 
signals are oceanic in origin. It is important to 
understand that these types of passive measurements 
permit generalization of accurate in-situ process 
observations to larger scales, but do not at this time 
have the accuracy of in-situ observations. Satellite 
oceanography currently is a highly descriptive en- 
deavor, only in limited cases will it be possible to ob- 
tain very quantitative results. However even in its 
presently crude quantitative form it does and will 
enable workers to test hypotheses of ocean circula- 
tion and transient response to forcing. 
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shelf waters, 0100Z, 4 April 1975. At the latitude of Lake Okeechobee, a nearly cloud- free at- : 
mospheric front occurs nen cooler, dry ee air to the north and warm, moist ¢ 


Introduction 


he great difficulty of acquiring requisite 

observational data plagues all efforts at ocean 
analysis and forecasting. In recent decades ocean- 
ographers and acousticians have acquired measure- 
ments that indicate kinetic energy at mesoscales (< 
24 hours and tens to hundreds of kilometers) often 
exceeds that of the synoptic mean motion by as much 
as two orders of magnitude.' The importance of 
mesoscale processes and the means to observe and 
model them now commands the attention of much 
modern ocean research. Our understanding of meso- 
scale processes in coastal or oceanic boundary 
regions has at present inadequate observational and 
theoretical bases. Oceanic processes at those scales 
generally defy prediction or even detection. Satellite 
scanning radiometers have the potential to supply the 
needed synoptic data but for two serious problems. 
First, the atmosphere, even without clouds, can 
obliterate the desired oceanographic information. 
Secondly, the radiometric temperatures may repre- 
sent only the skin temperature of the sea. We have 
found, and describe here, those times and places 
where these problems are minimized and a 
‘‘window’’ on the mesoscale structure of the upper 
ocean occurs. 

Over the continents in winter, the vast high- 
pressure systems of cold, dry, and dense polar air are 
generated largely by radiative cooling. The mobility 
of weather systems in mid-latitudes results in periodic 
outbreaks of this air across coastal regions into the 
adjacent ocean basins. The ‘‘window’”’ results from 
two aspects of this process. First is the development 
of clear, dry, relatively uniform air over the sea, 
through which best quality radiometric measure- 
ments can be made. Second is the deeply mixed con- 
ditions caused by strong winds and negative oceanic 
heat fluxes, which make surface measurements repre- 
sentative of much of the upper water column. Under 
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these conditions infrared and visible sensors of satel- 
lites can detect the mesoscale patterns of surface tem- 
perature and turbidity particularly well. These out- 
breaks occur in 3-6-day cycles, bringing one or more 
days of ideal imaging conditions at a time and yield- 
ing imagery of excellent quality. 


Atmospheric Effects 


Implicit in the above discussion is the fact that 
polar continental air is superior to maritime air for 
satellite observations and measurements. Satellite 
scanning radiometers collect the planetary radiances 
that reach altitudes of ‘‘inner space’’ (850 km for 
NOAA satellites). The radiances in the thermal in- 
frared part of the spectrum are a combination of 
radiant energy emitted from the earth, cloud, or sea 
surface that has been altered more or less by the 
intervening atmosphere. This alteration by a cloud- 
free atmosphere is a complex function of molecular 
absorption and reemission processes. In the simplest 
terms, a variable fraction of the surface-emitted in- 
frared radiation is absorbed and the energy is re- 
emitted at a different temperature, mainly by the 
triatomic.molecules of the atmosphere (HO, O;, and 
CO.). Of these, water vapor is the most effective and 
variable in distribution. Its effect varies along the 
radiation path as a function of its concentration and 
temperature. Under warm and moist conditions the 
sea-surface temperature patterns may be completely 
obliterated even without condensation and cloud for- 
mation. The entire absorption and reemission effect 
is, however, minimal in winter when cold, dry air 
overlies the lands and seas of mid- and high-latitude 
regions. 

The whole problem of attenuation of infrared 
measurements has been approached by most investi- 
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gators from the standpoint of correcting satellite- 
derived effective radiation temperatures to closely 
approach true surface temperatures measured dir- 
ectly by ships. Through all types of atmospheres, in- 
vestigators have found satellite temperatures lower 
by 1-10°C. A uniform cooling of the field of sea- 
surface temperatures is not so serious a problem for 
sea-going users.” The satellite data can be calibrated 
with surface measurements transmitted from ships or 
buoys to make full use of the finely resolved horizon- 
tal patterns of temperature. Cooling of radiation 
temperatures is, however, not the whole attenuation 
problem. Our studies of the satellite- and aircraft- 
based surface temperature measurements over 
coastal regions indicate that the atmospheric effects 
can be usefully separated into two problems: 


(1) Atmospheric attenuation—the lowering of 

radiation temperatures from the surface- 
measured values due to the imperfect trans- 
parency of the atmosphere. 
Suppression of temperature differences—the 
reduction of measured surface temperature 
differences between warm and colder parts 
of the scene caused by the wavelength varia- 
tion of radiative transfer through the at- 
mosphere. 


This second problem, from a practical standpoint, is 
more serious, since it means that atmospheric effects 
could obliterate oceanic fronts and mesoscale fea- 
tures of tactical interest, to mislead the user. 

A number of studies of atmospheric effects on 
infrared radiation measurements of sea-surface 
temperature have been conducted at aircraft alti- 
tudes. In June 1971 off the Massachusetts coast an 
aircraft radiometer was unable to detect temperature 
changes of 2°C in the sea surface with a moist layer 
of about 300 m depth present over the sea.* Another 
well-documented case was reported by scientists of 
the Research Division, Japan Research Committee of 
Environmental Remote Sensing from aircraft 
multispectral scanner flights over Ofunato Bay, Hon- 
shu, Japan. In the morning on 25 July 1978 the in- 
frared channel showed excellent-quality imagery, 
with absolute temperatures slightly cooler than 
ground truth, as expected. That afternoon, thermal 
contrast was almost nil, anc aircraft measurements 
of surface temperatures were warmer than surface 
truth measurements. A review of the experimental 
results revealed that the morning flights had been 
conducted after the nighttime land breezes had 
brought drier, cooler air out over the coastal bay. In 
the afternoon the sea breeze had brought warm, 
moist air in from the Pacific, which obliterated the 


gradients in the same manner as we have seen in 
satellite data. Radio-sonde measurements available 
simultaneously with the flights show that the 2 km- 
thick atmospheric column beneath the aircraft of the 
morning flight was an average of 2°C cooler and had 
1 g cm~? less precipitable water than the afternoon 
air mass. 

Atmospheric suppression of surface temperature 
differences is well illustrated with the NOAA-4 in- 
frared image (10.5-12.5 wm band) acquired at 0100Z 
on 4 April 1976 over the Florida Peninsula (Figure 1). 
The temperature gradients associated with the Gulf 
Stream are seen to disappear near the bottom of the 
image. The Florida Peninsula at this time was partly 
overrun by an atmospheric cold front moving rapidly 
south. The air mass boundary between the continen- 
tal air and the tropical maritime air, visible just south 
of Lake Okeechobee, was unusually cloud free. Sur- 
face observations indicate zero cloud cover for all 
stations, except for two-tenths registered at West 
Palm Beach on the east coast, which was directly in 
line with the clouds seen seaward on the image. This 
cloud mass is the landward extension of the cloud 
shield of the cold front visible over the Atlantic. Air 
temperatures are 20-21 °C with dew point readings of 
2-4°C just north of the front, and 26°C with dew 
point readings of 18-20°C to the south. There is a 
35% humidity drop and a 10 g m~3 moisture content 
difference in the air between the surface stations 
across the front (Vero and West Palm Beach 
stations). Analysis of the Miami radiosonde data 
before and after the frontal passages indicates that 
the warm, moist atmospheric conditions, with 2.6-3 
gcm~? of total precipitable water (TPW), were 
replaced by colder, dry air with 1.1-1.5 gcm~? TPW. 
Nearly all surface temperature differences at the edge 
of the Gulf Stream and Florida Peninsula have been 
suppressed in the radiation temperature field by these 
atmospheric conditions. 

Maul and Sidran‘ calculated that an 8°C 
temperature difference at the sea surface, when 
viewed vertically through the 30°N July Standard At- 
mosphere, would be recorded as less than 3°C by 
satellite and further reduced to 1°C when viewed 
through the atmosphere at an angle of 60° off the 
vertical. Using the Maul and Sidran model, the radi- 
ative transfer process through the two air masses over 
Florida at the time of the image (Figure 1) has been 
simulated, and the results are shown in Figure 2. The 
solid line is the distribution of atmospheric water 
vapor in terms of absolute humidity versus altitude. 
The dashed lines are plots of the calculated radiation 
temperatures, as they would be measured by a down- 
ward-looking radiometer, rising through the atrno- 
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sphere for surface radiation temperatures of 20 and 
22°C. Starting with the 2° temperature difference at 
the sea surface (20 and 22°C), both atmospheres 
reduced the temperatures and the temperature dif- 
ferences calculated for the top of the atmosphere. 
The warm, humid air reduced the temperature dif- 
ference from 2 to 0.8 °C, whereas the cooler, dry air 
reduced the temperature difference from 2 to 1.4°C. 
Examination of the moist air plot shows that the sup- 
pression of the temperature difference has occurred 
in the lower moist boundary layer of the atmosphere 
(lower 3 km). The calculated radiation temperature 
has been increased in the lower 3 km prior to the 
decrease in the upper atmosphere. This is a result of 
the infrared radiance contribution made to the over- 
all upward radiation by the warm, humid maritime 
air. The suppression of temperature differences is the 
cumulative effect of two processes. The Lambert- 
Beer Law states that the decrease in spectral intensity 
of radiation by absorption with distance is a function 
of the intensity of the initial radiation as well as of 
the concentration of the absorbers. The second proc- 
ess is the emission of infrared energy by the atmo- 
sphere itself. The result of both processes is to reduce 
the effective brightness temperature of a warmer 
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body more than that of a cooler body. 

Our field data suggest that atmospheric attenua- 
tion occurs under all conditions but that the suppres- 
sion of temperature differences is an air-mass- 
dependent phenomenon. We have attempted to 
quantify this effect by obtaining surface temperature 
measurements of two distinct bodies of water simul- 
taneously with an offshore atmospheric sounding 
and a satellite overpass. We accomplished this five 
times with the NOAA-S satellite in February 1977 in 
the coastal waters of the north-eastern Gulf of Mex- 
ico.* In the cases of higher levels of TPW, the ther- 
mal contrast between the warm and cold waters was 
suppressed, and the suppression appears to increase 
with increasing TPW. In the cases of cold, dry air 

over the sea, the atmospheric effect appears to be 
primarily one of uniform lowering of surface 
temperatures, without loss of gradients in the data 
field. The few observations available with the older 
lower resolution satellite radiometer provide only a 
preliminary look at the ‘‘water-vapor versus 
infrared-imagery contrast’ relationship which re- 
quires further study. A number of atmospheric 
radiative transfer models are available**, but we 
have found that they seriously overestimate the at- 
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Figure 2. Atmospheric humidity profiles and simulated radiative transfer of infrared energy for 20 and 
22°C water surface through the two air masses, 3 and 4 April 1975. 
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mospheric attenuation. They do-somewhat better in 
estimating temperature difference suppression. 


Some Mesoscale Structures of 
the Upper Ocean 


Examples of the best quality infrared imagery 
available from the new generation of NOAA 
satellites are drawn from the study of the Korea 
Strait and adjacent seas sponsored by the ONR 
Coastal Sciences Program. Perhaps the most severe 
outbreaks of polar continental air in the world occur 
over this oceanic region off northeast Asia. A great 
variety of marine environments occur here, ranging 
from the Kuroshio-dominated edge of the Pacific 
Basin to the shallow, runoff-freshened reaches of the 
Yellow Sea and East China Sea. The Kuroshio- 
Tsushima Current System complexly branches and 
distributes warm, saline waters to the North Pacific, 
Sea of Japan, and Yellow Sea after transport through 
and mixing within the East China Sea. The latter two 
basins have vast areas of shallow water of <100m 
depth, strong tidal activity, severe atmospheric forc- 
ing events, and among the largest freshwater and 
fine-grained-sediment -fluxes on the globe. In these 
shallow-water regions and adjacent shelf edges, hori- 
zontal mixing is extremely important as runoff from 
coastal regions encounters and is mixed with the 
equatorial waters of the Kuroshio-Tsushima system. 
The vigorous winter-season cooling processes 
destratify the water column, forming oceanic fronts 
between shallow and deep waters and waters of 
northern versus southern origin. These effects 
facilitate detection and tracking of various water 
types by surface temperature using satellite infrared 
scanners and by turbidity (color) using satellite 
visual-range sensors. 

Geographic orientation is provided by Figure 3, 
illustrating ocean basins, land masses, and national 
boundaries. Figure 4 is a daytime infrared image of 
the region from the older NOAA-S satellite of 8 April 
1978. It shows the patterns of earth-emitted radia- 
tion, with highest intensities (dark) over the Kuroshio 
and sun-warmed land masses, lowest intensities 
(white) over the cold cloud masses, and various gray 
tones showing sea-surface temperatures with well- 
developed mesoscale convolutions. The data from 
this older satellite system had a spatial resolution of 
0.9 km (at nadir), a 0.5 °C quantization interval, and 
an analog data link with a resulting noise level of 1-2 ° 
RMS. The current generation of NOAA satellites 
(TIROS-N series) has a spatial resolution of 1.1 km, 
a relatively noise-free digital link, and a 0.125°C 


quantization interval. The data are processed here to 
a 0.25°C resolution. Both systems have sun- 
synchronous orbits with a 12-hour repeat period. 

Data we have processed from one of the new 
generation of satellites (NOAA-6) over the Korea 
Strait are shown in Figures 5-9. These are parts of 
two overpasses from 8 and 17 November 1980 ac- 
quired in computer-compatible tape format from Na- 
tional Climatic Center and analyzed at the Remote 
Sensing and Image Processing Laboratory of Loui- 
siana State University. Data from the arrays of 
measured temperatures acquired by the satellite are 
presented in three formats: imagery, temperature- 
distance transects, and spot samples. Spot samples 
are the mean values of temperature within the tri- 
angular areas labeled on the imagery. The images 
were individually enhanced to maximize the gray- 
shade contrast between the oceanic features. Imagery 
in Figures 5, 7, and 8 were specially processed with 
very cold land and cloud areas electronically black- 
ened to improve photographic reproduction. Aside 
from these artificially blackened areas, the darker 
tones represent warmer temperatures. The white lines 
in Figures 5 and 6 replicate ship track lines from the 
Korea Strait-East China Sea-Yellow Sea Remote 
Sensing Experiment (KEYRSEX 80), which was con- 
ducted in November and December 1980. During this 
oceanographic experiment soundings of the ocean 
and the atmosphere were obtained for correlation 
with satellite imagery. The numbered transects in 
Figure 8 do not indicate ship tracks, but are pre- 
sented to show surface temperatures and the strength 
of the gradients observed in the image. 

Overviews of both scenes in Figures 5 and 6 
show the warm waters of the Tsushima Current ex- 
tending from the south to a strong front between 
Cheju and Tsushima Islands. During this experiment 
and throughout the year this front is convoluted into 
a series of mesoscale loops and meanders. Near 
Tsushima Islands the warm waters split into two 
branches. The main branch, the Tsushima Current 
proper, flows eastward along the coast of Honshu, 
Japan, toward the Tsugaru Strait. The smaller north- 
ward branch is the East Korean Warm Current, fol- 
lowing a broadly meandering path parallel to the 
peninsula. Radiation temperature differences across 
the front on both days are shown in the transects of 
line E for each. Sea-surface measurements and STD 
stations acquired along line E on 17 November 
verified the presence and strength of the front. Table 
I summarizes a comparison of ship and satellite mea- 
surements. 17 November was a relatively warm and 
humid day for this region at this time of year, with 
surface air temperatures of around 16°C and total 
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Figure 3. Geography of lands and seas near Korea Strait. 
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Figure 4. NOAA-S satellite thermal infrared image (10.5-12.5 mm) covering the Sea of Japan, Yellow 
Sea, East China Sea, and a small part of the North Pacific, 8 April 1978. A surge in the winter 
monsoon has brought cloud-free polar continental air out over the entire region. Clouds 

(white) occur in a broken mass over southern Japan and in a broad band to the southeast. 


precipitable water levels of 1.5 g cm~?. As noted in values, but the temperature of the warmest waters 
Table I, radiation temperatures of both water masses _ was affected the most. The lowest temperatures and 
are lowered by several degrees from actual surface salinities were found by surface measurements to oc- 


Table I 
| Coastal Region Tsushima Current Cross Front 
_ Waters Waters AT°C 


Ship 14°C (32.8 °/,.) 22° (34.3 °/..) 8°C 
Satellite 11.5°C 18.25°C 6.75°C 
Ship minus 
Satellite 2.5°C 3.75°C 1.25°C 
(The AT suppression) 
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Figure 5. 


NOAA-6 satellite thermal infrared image (channel 4, 10.5-11.5mm) of the south coast of 


Korea and Korea Strait, 2325Z, 8 November 1980. Mean temperature of the contained picture 
elements is labeled next to the triangular sample areas. Lettered lines are locations of ship 
tracks of cruise that commenced on 14 November. Line E of satellite temperatures versus 
distance is plotted below. Cold land and cloud tops have been electronically blackened to op- 


timize the display. 


cur just north of the front rather than in the coastal 
regions. This narrow belt of cool, low-salinity water 
extends from the East China Sea-Yellow Sea region 
west of Cheju eastward through the strait. It is shelf 
water known to discharge eastward between the 
north wall of the Tsushima Current and the local 
coastal water from late summer to early winter. 
Electronically expanded sections of these images 
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are presented in Figures 7 and 8. The narrows of the 
strait by the Tsushima Islands from 8 November are 
enlarged to show mesoscale details with greater clar- 
ity. The radiation temperature patterns and values 
from the southern Sea of Japan, just east of Korea, 
from 17 November are shown in Figures 8 and 9. The 
extreme enhancement of the image provides water 
mass outlines and temperature stream line patterns. 
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The land masses are extremely cold at this hour (7:30 
a.m. local sun time) and so produce an abrupt tem- 
perature drop at the land-water boundary. Low 
cumulus clouds to the right of the imagery near 
Japan generate the distinct sawtoothed pattern seen 
in transect 1 of Figure 9. From this perspective, 
Tsushima Island mass acts as a fixed water vane, 
deflecting parcels of Tsushima Warm Current waters 
from the northward meanders of the frontal bound- 
ary along the path of the East Korean Warm Cur- 


Figure 6. 


rent. As a result, alternating warm and cold water 
parcels flow north along the east coast of Korea. 


Conclusions 


A practical consequence of these studies is the 
knowledge that usable information on the horizontal 
temperature structure of the upper ocean may be rea- 
sonably expected behind outbreaks of polar conti- 
nental air. Under other conditions, however, such as 


An electronic expansion of the vicinity of Korea Strait from the previous 8 November overpass. This 


subscene includes Korean coast to teh northwest, Tsushima Islands in the center, and clouds plus other 


islands to the south 


t. Sea-surface temperature patterns show Tsushima Warm Current 


waters interleaved with the cold waters of the coastal region. 


NR Reviews 


L é 
> | 
— 


AA-6. CH-4 
W270. 14 
@STAL STUDIES IN 
JECTRICAL ENG. 
MUIS IANA STATE 
WR CODE 462 


THE 


NHY 


DEPT. 
UNIV. 


Ww 

¢ 

rf 

= 

= 

= 


12.00 


10.00 


80.00 


in humid tropic regions or deep ocean areas where 
maritime air masses predominate, the single-channel 
data from a satellite infrared radiometer (8-13 um 
band of the spectrum) has little chance to faithfully 
record surface temperature patterns. On longer time 
and space scales, the atmospheric marine boundary 
layer (MBL) is strongly influenced by the sea surface, 
its major source of heat and water vapor. Satellite 
observations of radiation temperature patterns may 
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Figure 7. NOAA-6 satellite thermal infrared image (channel 4, 10.5-11.5 tm) of the Korea Strait and 
parts of adjacent seas, 2326Z, 17 November 1980. The temperature-distance plot corresponds 
to a line of oceanographic measurements made by a ship on the same date. 
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follow oceanic fronts in humid regions by virtue of 
the strong influence that surface conditions have on 
the MBL. As sea-surface features move or change, its 
thermal counterpart in the MBL may lag behind, or 
vice versa. We would judge that cloud-free temper- 
ature patterns in warm, humid oceanic regions would 
not serve as an operationally reliable link to the up- 
per ocean because of the strong surface atmospheric 
absorption and emission effects noted above. 
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Figure 8. NOAA-6 satellite infrared image of the Korea Strait and southern Sea of Japan electronically 
expanded from the 17 November 1980 overpass. Flow paths of the Tsushima and East Korean 
Warm Currents are shown in the sea-surface temperature patterns of the southern Sea of 
Japan. Temperature transects 1 and 2 are presented in Figure 9. 


However, analogous to cold-air outbreaks are 
localized mesoscale winds such as the mistral and 
land breezes. Although they are smaller in spatial and 
temporal scales, they too blow from land to sea. The 
drier and cooler air characteristic of such offshore 
winds may also provide ‘‘windows’’ for remote sens- 
ing of coastal and continental shelf waters. Certainly 
more meteorological/oceanographic/remote sensing 
team studies are needed to explore and make prac- 
tical the use of atmospheric ‘‘windows”’ on the sea. 

Avoidance of the atmospheric limiting effects is 
perhaps possible using other parts of the spectrum or 
using a multispectral approach, with calculations 
from spectral attenuation differences to correct each 
image element. The mid-infrared data from channel 
3 (3.5-3.9 wm) of the NOAA satellites is much less 
susceptible to atmospheric effects, and preliminary 
studies indicate promise for this sensor.*’’ It is usable 
at night only, because of the strong solar component 
in this part of the spectrum. At this time it appears 
that true all-weather, day-night sea-surface temper- 


ature data must await development of high-resolu- 
tion multispectral passive microwave radiometers. In 
the meantime there is much oceanography to be 
learned, models to be tested, and operational support 
possible under the outbreaks of polar continental air 
around the globe. 
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Figure 9. NOAA-6 temperature-distance transects 1 and 2 from Figure 8, 17 November 1980. The north 
wall of the Tsushima Current near Japan (upper) and the west wall of the East Korean Warm 
Current (lower) are 2.5°C radiation temperature steps. Intersection with land and cloud tops 
produces the lowest temperature values. 
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Field Study of 


“Anomalous Gray-Shades”’ 


S atellites have been used for many years to 
observe clouds and cloud systems to better 
understand atmospheric processes. Clouds are 
rather easy to ‘‘see’’ with a satellite due to their 
high reflectance or albedo. Recently introduced 
satellite sensors are capable of detecting variations 
of very low albedo in the atmosphere. These 
variations of brightness away from bright cloudy 
regions are called ‘‘anomalous gray-shades’’ and 
are frequently observed in visible and _near- 
infrared images from other satellite systems. 

A team of investigators under the direction of 
Professor Thomas Vonder Haar at Colorado State 
University, funded by the Office of Naval 
Research’s Selected Research Opportunity 
Program has been studying these anomalous gray- 
shades. 

To investigate these features, a coordinated 
research effort was conducted between 19 Sep- 
tember and 2 October 1980, off shore of Southern 
California. Airborne state-variable, air motion, 
particulate and radiometric measurements were 
obtained to determine the structure of gray-shade 
regions identified in satellite images. Preliminary 
results indicate that certain of these features are 
associated with the marine boundary layer. They 
seem to be related to variations in size and con- 
centrations of large particles primarily within the 
boundary layer. (Figure 1 shows a set of features 
observed with the NOAA-6 Advanced Very High 
Resolution Radiometer (AVHRR). They were 
observed on the 28th of September in the wake of 


Anomalous Gray-Shade Features in the Wake of 
Santa Rosa and Santa Cruz Islands off of Santa 
Barbara, California as Observed by the NOAA-6 
Advanced Very High Resolution . Radiometer 
(A VHRR) on 28 September 1980. 

Santa Rosa and Santa Cruz Islands just off of 
Santa Barbara, California. The features are 
aligned with the boundary layer wind and 
although some of the brightness variations can be 
attributed to surface roughness variations, most 
of the signal is due to atmospheric reflection.) 
Anomalous gray-shade regions, when fully un- 
derstood, may provide great observational insight 
into stratocumulus formation and dissipation, 
humidity variations, visibility determination and 
forecasting as well as boundary layer dynamics 
and air-sea interactions. 


Figure 1. 
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mages formed from radar sea surface scatter 

return to aircraft and satellite platforms show 
coherent patterns that can be associated with current 
shear, eddies, internal waves, and patterns on the sea 
surface and bottom topography. Seasat SAR (Syn- 
thetic Aperture Radar) data show a wealth of de- 
tailed patterns. An example of this in Figure 1, shows 


a digitally-processed image of Nantucket Island, 
Massachusetts and its surroundings. The magnified 
part in Figure 2 indicates a velocity shear at the edge 
of a current of less than one knot running to the 
WNW, south of the Shoals. One cannot explain these 
changes in the wave packet structure of sea surface 
waves in terms of the results of the customary ap- 
proximation in terms of waves of infinitesimal ampli- 
tude or linear theory. The patterns in wave amplitude 
and wavelength are due to effects of finite wave am- 
plitude, where nonlinear phenomena become impor- 
tant. The tendency for waves to undulate in ampli- 
tude and frequency is typical of finite amplitude 
waves. We refer to these as amplitude (AM) and fre- 
quency (FM) modulations, and, since they are due to 
nonlinear dynamics, we call them nonlinear wave 
modulations. 

I shall attempt to describe in this article some 
features of modulational dynamics for surface 
waves, and how the detection process used for SAR 
and Side Looking Airborne Radar data rely on inter- 
modulations of the scattering wavelets and the under- 


Surface Wave Modulation 
Patterns in Radar Images 


by 


Erik Mollo-Christensen 
Massachusetts Institute of Technology 


lying dominant wave field. Then I will go on to sug- 
gest that the same kind of forced modulation and in- 
termodulation changes the patterns in the dominant 
wavefield to make eddies, shear layers and other 
coherent motions visible. First consider modulations 
of the surface wavefield. 


Surface Wave Modulations 


Both surfers and boatmen know that swell and 
wind waves vary in amplitude and arrive in packets 
or groups. A wavefield tends to develop variations in 
wave amplitude, a, because, at a given wave steep- 
ness, a dip in amplitude (a,, > 0) tends to move 
faster than a positive bump in wave amplitude (a,, < 
0). So, dips and bumps will tend to separate rather 
than merge and cancel. This instability process makes 
disturbances in wave amplitude grow initially, when 
a uniform amplitude wavefield is disturbed. But near 
the bumps, there tends to be larger wave steepness 
ak = 27a/h, defined in terms of amplitude, a, and 
wavelength, A. Steeper waves tend to propagate 
faster, so there will be a tendency for the front of 
wave groups to steepen further. 

These effects of amplitude variation and steep- 
ness on the phase speed of finite amplitude periodic 
waves (Stokes waves) are contained in the dispersion 
relation, relating phase speed c, wave number k = 
2m/d, and wave steepness ak = € 
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= gk[] + + a,,/(8ek)]. 


Eventually, in the interplay between steepness and 
bumpiness effects, an equilibrium situation is estab- 
lished, so the mature wave field consists of constantly 
changing groups of waves, maintaining an average 
steepness and wave amplitude. This is the typical sea 
surface wave field under weakly generating wind 
wave conditions, with nearly constant average prop- 
erties, but always changing locally in terms of group 
lengths and amplitude. Now and then accidents 
happen, in the sense that a further instability occurs 
where subharmonic transitions take place, and a crest 
is lost. Due to such accidents, the average wavelength 
of a wavefield will tend to increase and the wave 
steepness decrease until there is an open ocean swell 
wavefield of modest steepness. The description I 
have given is based on analytical results for nonlinear 
waves, but not all the detailed steps have been veri- 
fied as occurring in the order described, so, in that 
sense, my description is somewhat speculative. 


Figure 1. Digitally processed image of Nantucket Island and the 
surrounding waters obtained from SEASAT SAR data. 


Linear and Nonlinear Wave Dynamics 


Linear theory is excellent for describing wave 
development over limited distances and times; for 
surface waves, linear approximations are good for 
propagation distances, L, of much less than )/e’, for 
larger distances nonlinear effects become more im- 
portant. For 100 m wavelength swell of amplitude 1 5 
m or nearly 5 ft. crest to trough height, « = 0.1 and 
the propagation distance is 10 km. For small wavelets 
of 10 cm wavelength and 0.5 cm amplitude e = 0.3, 
the propagation distance is only one meter. 

For propagation over longer distances, the cu- 
mulative effects of nonlinear terms in the equations 
of motion become important, and one has to develop 
the dynamics to third order in wave steepness and ob- 
tain the Davey-Stewartson equation' in wave- 
complex amplitude A(x,y,t), 


i(A, + c,A,) 
Fp: (Ag — 2Ay) (2) 
- SARA = 0. 


Figure 2. Magnified portion of Figure 1., showing the shear 
instability of a current edge shear layer. 
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The water surface deflection is 


z = $(x,y,t) = RefA(x,y,t) (3) 


Here, k and w are the average wavenumber k = 27/\ 
and average frequency w = 22/T, T is the wave 
period. 

The Davey-Stewartson equation gives a valid ap- 
proximation over distances of propagation much 
smaller than L; = )/e?*. (L; is not the same as the L 
for linear approximation. It gives a ‘‘good’’ descrip- 
tion for propagation time scales of order T/e? and for 
group lengths of order \/e. For the swell example, 
h = 100m, « = 0.1, L; = 100 km. For the wavelets 
of \ = 10cm ande = 0.3, L; = 30m. 

The Davey-Stewartson equation (2) has exact 
solutions? in the form of periodic and solitary distur- 
bances. In two dimensions, this equation reduces to 
the cubic Schrodinger equation’, first used in wave 
mechanics to solve for the electron as a wave packet. 
Also when applied to sea-surface wave modulations, 
the solutions to the cubic Schrodinger equation have 
quantum properties, limiting the energy content of a 
single wave packet for a given wavelength. The solu- 
tions also possess certain conservation properties, so 
that wavepackets that collide will survive collisions. 
This suggests that wave packets can have long half 
life, and that wavepacket descriptions of a sea- 
surface wave field may be useful for some purposes.* 
But, as will be argued further along, the modulation 
field is also very sensitive to disturbances, which can 
alter the wave-packet structure of a wave field en- 
countering current, wind or topography. 


Forced Modulation 


The modulation field, as described by equation 
(2), will, when otherwise undisturbed, reach an 
equilibrium state where the linear and nonlinear 
terms balance. A current shear dv/dx = vk, where 
the mean velocity changes by an amount v in a wave- 
length d, will introduce additional terms in the 
Davey-Stewartson equation. One term comes from 
the energy-flux divergence term in the energy equa- 
tion’, and introduces a term (vA), ~ (vAk) in equa- 
tion (2). There are additional terms due to radiation 
stresses’, but they are of the same order-of- 
magnitude. 

Velocity shear will be as important as the non- 
linear term in equation (2) when 


Avk = 


2 
2 
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setting 


c, = = and k? |A|? = e’, 


one finds that 


Vv 
is the criterion for current shear to affect modula- 
tions. This shows that small current shear can be im- 
portant. 

For an example of a 100 m wavelength swell of 
steepness « = 0.1, the group velocity is c, = 6 m/s 
and the current variation required to disturb the 
modulation field is v = e*c, = 6cm/sec = .12 knots. 
For the 6 cm waves of steepness « = 0.3, the current- 
change needed is of order v = 1.5 cm/sec. For less 
steepness, with e = 0.1 and X = 6 cm, we find the 
current-variation needed to upset the modulation 
field to be of order v = 1.5 mm/s. 

This demonstrates how very small coherent dis- 
turbances in the current field can cause disturbances 
in wave modulation patterns in wavefields of moder- 
ate steepness, and if the disturbances are coherent, 
the modulational disturbances will also be coherent. 
The steeper the waves, the larger the disturbance 
needed to affect the modulation field. But ship wakes 
and weak eddies will show up in the modulation field 
of short wavelets of \ = 6 cm under a wide range of 
velocity disturbances down to a few cm/s. 

Disturbances in wave propagation speed due to 
bottom topography will affect the modulation field 
the same way; for shoaling from very deep to a depth 
of 20 m will change the group velocity of 100 m 
waves from 6.25 m/sec to 5.32, so that the change is 
v = 0.93 m/sec, sufficient to change the modulation 
field. This suggests the reason why one can see the lo- 
cation of the continental shelf edge in some SAR im- 
ages, and why topography shows up so very clearly. 

So far, we have discussed how waves, as they 
propagate, eventually accumulate effects of nonlin- 
earity until they add up enough to change the wave- 
field by introducing modulation in phase and ampli- 
tude. The steeper the waves are, the quicker they 
show the effects of nonlinear processes. But through 
all this, the linear approximation is a good one for 
small time and distance of propagation. So observa- 
tions of wavefields over short distances and times will 
tend to be dominated by linear effects and reveal lit- 
tle of the nonlinear dynamics. But for the much 
larger scales of observation made possible by remote 
sensing, the nonlinearities have room to show up, 
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and give disturbances their chance to modify wave 
modulation fields. 


Intermodulation 


The current that disturbs the wave modulation 
field may be a fluctuating current due to other waves. 
A field of \ = 10 cm waves riding on waves of 100 m 
wavelength of 3 m amplitude (« = 0.10) will be sub- 
ject to current fluctuations of + 1.2 m/sec. Accord- 
ing to our calculations, this current is sufficient to 
dominate the modulation field of the 10 cm waves. 
The orbital velocity field is convergent-divergent, so 
that it plays accordion with the wavelet field, induc- 
ing large steepness and leading to fourth-order ef- 
fects becoming important, the wavelet field may 
undergo subharmonic transitions and even wave- 
breaking at certain phases of the underlying domi- 
nant sea. 

The mechanism of forced modulation of the 
wavelets by the dominant sea is the basis for the SAR 
detection of sea surface waves**:'°. The SAR signal 
is scattered from wavelets and the intensity and Dop- 
pler shift of the reflected wave is interpreted as a 
modulation signal due to the dominant wave and, 
through a suitable algorithm, translated as sea sur- 
face shape due to the underlying waves. This inter- 
pretation algorithm for SAR and other radar scatter- 
ing imagery is, of course, a crucial element in data 
interpretation. Because of the complexity of the in- 
termodulation processes, one has to develop better 
algorithms by simultaneous observation of both 
wavelets and underlying waves over a range of wind 
and wave conditions. This should preferably be done 
without the need for flying radars on aircraft or 
satellites, since the cost of sufficiently extensive 
overflights would be prohibitive. I have therefore 
designed a scattering-wave gauge which simulates 
radar backscatter. 


The ‘‘Harp’’ Wave Gauge 


A student, J. Menoher, and I have built and 
tested a simple wave gauge whose output simulates 
radar scatter return. The wave gauge consists of a 
single capacitance wire, strung on a frame with 
parallel vertical segments a distance f apart, as shown 
in Figure 3. The capacitance of the wire varies pro- 
portionally to the sum of the water level displace- 
ments at all the wires. A capacitance detector circuit 
will then provide a wave gauge output proportional 
to the sum of displacements. For very long wave- 
lengths, compared to the length of the array, the 
gauge will work like a single point wave gauge. For 
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“‘Harp”’ wave gauge, with a wire strung in parallel runs. 
Resonant wave contour shown. 


Figure 3. 


waves of wavelength equal to the wire spacing, the 
capacitance variation of all the wires will be in phase, 
and the output will be a maximum. For wavelengths 
twice the wire spacing, the output will be zero. The 
output amplitude will vary with wavelength as shown 
in Figure 4. The wave gauge functions as a resonant 
scattering detector. It can also be built for less than 
the cost of a satellite; its capital cost will be com- 
parable to an aircraft flight hour. 

A first test of a ten-wire wave ‘‘harp’’ gave the 
results shown in Figure 5. The upper trace is the 
‘tharp’’ output and the lower one is the output of a 
single-wire capacitance wave gauge placed next to it 
in a wave tank. The frequency of the wave-generating 
paddle increases slowly, and the resonance peak is 
evident. But it is also interesting to note that there are 
bursts of high-frequency output for incoming waves 
of frequencies of 1/2 and 1/3 of the frequency of the 
resonant wave. 

This means that a scattering detector also detects 
higher harmonics of Stokes waves. The same must be 
the case for a Bragg scatter-return of radar signals. 
This result illustrates the complexity of interpretation 
of scatter returns, and indicates the need for careful 
observations to establish relationships between scat- 
tering waves, wind and dominant sea. Data process- 
ing of data from such an experiment or from SAR 
data, for that matter, requires the use of processing 
methods that produce information about modula- 
tions in phase and amplitude and intermodulations. 


Eduction of Modulations from Signals 
What we need is information about the complex 


modulation variable A in equation (3). One method, 
taken from communications engineering, uses 


Hilbert Transforms'' and has been used for surface 
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Figure 4. Calculated output of ‘‘harp”’ scattering wave gauge a function of wave length. 
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Figure 5. Output of “‘harp’’ wave gauge (upper trace) compared to signal from a single wave gauge 
(lower trace) as wave frequency is slowly increased showing higher harmonic resonance 
(A&B) and main resonance peak (C). (Result obtained by J. Menoher). 
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wave processing by my colleague W. K. Melville. One 
of our students, G. Sahar,'? has used the method to 
process laboratory wind wave data. 

Figure 6 shows a wave gauge output, with the 
amplitude |A| superimposed, and, in the lower 
trace, the phase e of A = Re'®; the linear trend due 
to the carrier wave exp[i(kw-wt)] has been taken out. 
Once one obtains A, one can relate it to the other 
scale processes that drive or interact with the modula- 
tions. The interpretation of SAR return by scattering 
wavelets and relating it to the dominant sea is a dif- 
ficult problem and is similar to the challenges of 
relating wave modulation to internal waves and other 
large scale processes. They present similar data proc- 
essing problems. 

In learning to interpret data in terms of such in- 
termodulational relationships, we need to do some 
controlled experiments, guided by existing work on 
nonlinear waves and intermodulation and with an 
awareness of the limitations of simple and convenient 
approximations. 


Applications to Oceanography 


The real advantage of active electromagnetic 
sensors in Oceanography is their insensitivity to 
clouds and fog. One does not realize the prevalence 
of cloud and fog patches before one starts looking 
for visible and infrared satellite data for the sea sur- 
face. For some locations, such as the Drake Passage, 
for example, one may find one day a year with a 
good view. Other areas are nearly as frustrating; 
climatological records of sunshine, cloud-cover and 
fog are readily available to demonstrate this. Also, 
interesting events, such as hurricanes, gales and in- 
tense air-sea interaction do not occur on sunny days. 
So, without sensors that sve through clouds, we can- 
not obtain the synoptic view of patterns that are 
uniquely available from remote sensing devices. 

Good distance-resolution, another crucial fea- 
ture of SAR and SLAR data, derives from the fact 
that distance is measured in terms of time for a signal 
to return rather than through optical resolution of 
geometry. A spatial resolution better than fifty 
meters makes it possible to see the dominant surface 
waves and swell, internal waves and current shear. 
Because the dominant waves can be resolved, it is 
also possible to educe from the data the amplitude 
and modulation field of those waves, and to interpret 
the modulation patterns in terms of larger scale 
phenomena. It is worth emphasizing that, because 
the modulation field is governed by nonlinear proc- 
esses, One cannot use linear filtering theory to 
calculate resolution. The resolution of an image 
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needed to discover a 100 km eddy from the surface 
wave modulation field must be sufficient to resolve 
individual waves. 

As an alternative, one can build a modulation 
eduction operator into the processing scheme. One 
reason for the clear view of amplitude modulations 
afforded by existing SAR imagery is that somehow 
such an eduction operator is present in the data proc- 
essing chain, possibly by inadvertence or because of 
the skewness of modulations of scatterers on surface 
waves. 

I will next give a few examples of data use in 


oceanography. 


Density Structure from Internal 
Wave Patterns 


When and if SAR and similar active sensor sea 
surface data become available, it will be possible to 
put into routine use a method for estimating vertical 
density structure that a student, A. da. S. Mascaren- 
has, Jr., and I developed a few years ago. We used 
Landsat images that showed patterns that could be 
identified with internal waves propagating on the 
seasonal thermocline. These internal wave packets 
occur in the New England continental shelf region. 
The waves are generated by tidal flow over topog- 
raphy, and the packets are generated once per M, 
tidal period, which is approximately 12.4 hours. One 
sees packets generated several tidal periods apart in 
many Landsat images. The spacing between packets 
is the distance travelled in 12.4 hours. This gives the 
propagation speed, c, of the wave packets. 

The linear approximation for internal waves on 
a density jump Ag/g at a depth h from the surface in 
otherwise deep water gives 


2= ph Oe. 
c? = gh 
Since we have found c from the satellite image, one 
can calculate h s If one knows Ag/e from other 
sources, one can find h. We argued that the density 
jump on the seasonal thermocline was due to the tem- 
perature difference, so that Age/e = — aAT. We 
then said that hAT was a measure of heat storage in 
the upper layer, and calculated its distribution for the 
months of July and August'*. (Figure 7). 

We were constrained by limited data availabil- 
ity, caused by fog and the Landsat data acquisition 
schedule. This method may well prove applicable to 
internal waves on the main thermocline, since one 
should expect such waves to be generate? all along 
the continental edges by the tides. We may need bet- 
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Figure 7. 


Heat storage in upper ocean derived from internal wave 
packet spacings. Mollo-Christensen (1980). 
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ter methods of processing radar data to see the sur- 
face modulation patterns generated by these internal 
waves. The distance between wave packets will give 
the speed, and one can again observe . Another 
phenomena that can be excited by the tides is shear- 
layer instability. 


Estimating Current Speed from 
Satellite Data 


Another example of interpretation of patterns 
visible in radar images is given by regarding the un- 
dulations of a current edge shear layer as a tidally 
forced hydrodynamic instability. '* 

Where a current passes a prominent topography 
feature and speeds up locally, oscillations and tur- 
bulence tend to decrease in comparison to the orga- 
nized motion; in part, this is due to the increase in 
mean velocity, without a corresponding increase in 
oscillations and turbulence. Also, since the current 
touches a boundary, there will be little oscillation 
normal to the boundary. Further downstream, where 
the current loses the support of the boundary and the 
current edge becomes a free shear layer, the flow 
becomes unstable. An unstable flow tends to amplify 
fluctuations due to external disturbances, so that the 
shear layer develops lateral fluctuations in the form 
of travelling waves of increasing amplitude. This is 
clearly visible in TIROS-N infrared images of the 
Gulf Stream near Cape Hatteras, North Carolina, 
shown in Figure 8. Another current edge can be seen 
in the Seasat SAR image of the Nantucket Shoals, 
where the edge of the WNW current south of the 
Shoals, reported in the U.S. Coastal Pilot, forms an 
undulating shear layer edge, shown in Figure 2. 

The excitation for shear flow instability will be 
velocity fluctuations caused by incoming internal 
waves, storms and tides. The most reliable excitation 
is the tide, and, in the examples shown, the M, tide. 

We can calculate the speed of the shear layer 
waves by measuring the wave length in the satellite 
image and dividing by the period of the M; tide. For 
the Gulf Stream just beyond Cape Hatteras we find 
c = 1 m/sec. Next, we use the fact that shear layer 
waves tend to move at a speed c = (U, + U,)/2, the 
average of the current velocity on the two sides. We 
now assume that the current velocity on one side of 
the shear layer is zero: U, = 0; one then finds 


2 (wavelength) 
(M-period) 


U, = 2C = 


For the Gulf Stream, we found U, =~ 2 m/sec, very 
close to the reported value. 

In order to verify such methods of current esti- 
mation, one needs a series of sea-surface images 
spaced a day apart or so, so that one can follow a 
shear layer wave. But regions where a cold or warm 
current passes a headland tend to be active mete- 
orologically, breeders of storms, regions of cumulus 
convection or just foggy. It is hard to find infrared 
radiometer images of Hatteras or the Nantucket 
Shoals free of clouds. The same goes for most other 
major currents like the Kuroshio, the Antarctic Cir- 
cumpolar, the East Greenland, and the Somali Cur- 
rent as well as the many minor currents we could 
analyze this way. 

Again, active sensor data, such as SAR or SLAR 
data, would be very important for estimating cur- 
rents, analyzing the growth of fluctuations and 
meanders as the first stage in formulating predictive 
models of such undulations. 


Ship Wakes as Diagnostic Tools 


Ship wakes show up in high resolution images 
such as those produced by SAR. The deformation of 
a ship wake by current shear can be used to deduce 
the shear, provided the ship keeps steady course and 
speed, and one knows the speed. One can find the 
speed from the wavelength of the stern wave. As 
anybody who has spent time leaning on the stern rail 
of a ship knows, the stern wave moves at the speed of 
the ship; this means that the stern wave propagates at 
the speed of the ship, v. If \ is the wavelength of the 
stern wave, one has the relation 


V = 


Combining the analysis of shear instability with use 
of ship wakes as a diagnostic can aid one in mapping 
out current fields and fluctuations in a region. This 
requires high-resolution sea-surface data without 
cloud or fog clutter. 


The Data Volume Problem 


High resolution data represent high data density 
per unit area of the sea surface. SAR data with a 
resolution of 15 x 15 m? contains approximately 30 
megabytes for a 100 km? data field. Transmission, 
storage and processing of such large data volumes are 
just now becoming routinely possible at manageable 
cost. The technologies of multiprocessors and mass 
storage devices are advancing so rapidly that within a 
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decade, data handling problems and costs will have 
changed for the better, and we must try to anticipate 
technological advance as best we can, and guess at 
where future bottlenecks will occur. 

Data archives and data distribution to research 
workers may end up on high-density disk media, 
while operational users will get their data via com- 
munications satellite. 


Improvement in Quantitative Interpretation 
of Data 


A crucial element in interpretation of sea- 
surface radar-scatter data are nonlinear wave proc- 
esses. We need to know such processes better, and we 
need to be able to analyze the data field to take apart 
the results of nonlinear interactions. We can learn 
much about nonlinear waves from laboratory and 
field experiments and much of the knowledge we 
need can be obtained without having to fly satellites 
or aircraft. 

The analysis methods we now have work with 
time series and joint time series, that is for functions 
of a single independent variable. We need to extend 
such methods to two-dimensional and three-dimen- 
sional data fields, such as those we obtain by remote 
sensing. 


Figure 8. 
AVHRR Data). 
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The nonlinear methods for time series are not 
always easy to generalize. One need only realize that 
while the order in which one applies linear processors 
to a signal may not affect the result, the order of ap- 
plication of nonlinear processors can make a large 
difference. For example 


sinJu(t) # Vsinu(t) andlogu(t) # (log u/t)?. 


Conclusion 


Radar scatter data are now useful in operations 
and research. Improvement in quantitative inter- 
pretation of the data requires research on wave 
dynamics and data processing. Technology of data 
transmission and storage will keep improving and 
will be available to handle remote-sensing data when 
needed. 

Patterns in remote-sensing images can be inter- 
preted quantitatively to produce information on cur- 
rent speed and density structure. For such use, the 
all-weather capability and the good spatial resolution 
of active sensor data will be especially valuable. Hf 


The edge of the Gulf Stream beyond Cape Hatteras, showing shear waves (from TIROS-N, 
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Introduction 


he evolution of computer technology has pro- 

vided scientists and engineers with the capability 
to process information efficiently in the form of im- 
ages or pictures. This has resulted in a great deal of 
research interest in image processing as exemplified 
by the May, 1981 special issue of the Proceedings of 
the IEEE' devoted to papers on image processing. 
The collection of papers emphasized two aspects of 
image processing: (1) the processing of image data to 
achieve desired improvements, and (2) image under- 
standing which is concerned with the way humans 
visualize and utilize image information. 

In the meteorological satellite community, scien- 
tists now routinely use computers to process image 
data to improve image quality and image understand- 
ing. For example, an algorithm to remove white 
background due to snow provides a better quality 
picture for cloud studies. Another example is the 
assignment of separate colors to cloud heights which 
allows the user to differentiate clouds quickly at dif- 
ferent levels and thus provides improved image under- 
standing. The two examples deal with flat images as 
viewed from the satellite. Currently, the majority of 
satellite imagery is viewed in this context and the user 
must visualize the three-dimensional structure of the 
atmosphere in his mind. 

Given a flat satellite image, contouring algo- 
rithms can be employed to operate on the visual or 
infrared data to generate three-dimensional cloud im- 
ages. Another alternative is to use polygons to ap- 
proximate a three-dimensional cloud surface. These 
methods result in three-dimensional surfaces that 
more realistically depict what a pilot will see from an 
airplane or what an observer will see from the 
ground. The result is improved image utilization. We 
will define image utilization as a transformation that 
results in a product that directly depicts a desired 
phenomenon. 
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This paper describes the use of satellite pictures 
and computer-driven imagery to represent clouds as 
three-dimensional objects in a viewer’s frame of ref- 
erence. Given a flat satellite picture, the procedure 
for generating the cloud pictures is described. The 
use of shading to provide a more effective three-di- 
mensional picture is illustrated. Then dimensioning 
of the images is described, again from the viewer’s 
frame of reference. 

The visual and infrared image data were col- 
lected by the ground station at the Atmospheric 
Science Department at Colorado State University. 
The computing and display equipment consisted of a 
Digital Equipment Corporation PDP-11/60 mini- 
computer and a RAMTEC 9050 Color Image 
System. The RAMTEC has a display of 512-by-512 
pixels with each pixel storing an 8-bit value. A video 
lookup table is used in the RAMTEC to define the 
color associated with ech possible pixel value. A 
13-bit word stores each color definition. The lower 12 
bits are divided into 3 groups of 4 bits each. Each 
group controls the intensity of one of the three stan- 
dard color guns in the monitor. This provides a 
possibility of 16 intensity levels for each color gun 
and a total of 4096 colors and intensities. Only 256 
colors and intensities may be defined at one time, 
however, and only 16 intensities can be generated for 
a single color. A more detailed description of the 
work is given by Johnson.? 


Two Dimensional Representations 


The data that were received from the satellite 
consisted of visible and infrared data. In the visual 
data, the pixel values are related to brightness which 
does not directly relate to altitude. This limits the 
usefulness of visible data for representing clouds in 
three dimensions. The infrared temperature data can 
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Figure 1. 


be correlated to altitude via data collected from bal- 
loons so that reasonably accurate altitude represen- 
tations can be achieved. 

Thus, the first task is to use balloon data to re- 
late each infrared pixel value to altitude. After this is 
done, a new image is generated to scale the altitude 
values between the values of 0 and 255 which repre- 
sent the lowest and highest altitudes respectively. The 
resulting new image can then be displayed again as a 
flat image. 

An example of the results is shown in Figures 1, 
2 and 3. Given the original image as shown in Figure 
1; a histogram plot of the data after conversion is 
shown in Figure 2. The relationship between pixel 
values and altitude is given by 


A= (64.5)I (1) 


where I is the pixel value and A is the altitude in 
meters. The resulting reconstructed flat image is 
shown in Figure 3 with a gray scale from 0 to 200. 
Two hundred was used as the upper bound since 
values above 200 are constant. A pseudocolored 
representation of Figure 3 is shown in Figure 4. This 
dramatically brings out the two clouds in the center 
of the picture and these are the two clouds we will 
place into a three-dimensional image in the next sec- 
tion. 


Original infrared image. 


To summarize, given a flat infrared image from 
the satellite, a conversion to altitude is made using 
the best available information to convert temperature 
values to altitude. Histogram values are then ad- 
justed between 0 and 255. This may cause distortion, 
but the impact of this is still not known. Then a new 
image is constructed and pseudo-colored if it is to 
show relevant features. 


_ 160 192 224. 256 
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Figure2. Pixel histogram plot; pixels converted to altitude. 
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Figure 3. 


Three-Dimensional Representations 


Given the infrared images, clouds in a desired 
sector can be converted to three-dimensional objects 
that can be viewed from any angle. The resulting im- 
ages can be shaded and dimensioned to provide the 
viewer with a realistic view of a cloud or group of 
clouds from any desired viewing angle. 

A surface can be created from the reconstructed 
data by means of a collection of surface drawing 
subroutines that are available from National Center 
for Atmospheric Research (NCAR)? in Boulder, Col- 
orado and were modified at CSU for use with our 
computer equipment. With the PDP-11/60 and 
RAMTEC system we are limited to an array of 
32-by-32 points. Since the number of pixels in the 
portion of the image used to create the surface is 
much greater, a method to utilize the data was 
needed. The method which we chose was to select 
pixels at regular intervals so that they form 32 rows 
of 32 columns within a selected portion of the image. 
The resulting surface can be generated as a true con- 
tour map or using an X, Y, Z coordinate system. 
Modification of the programs resulted in surfaces 
with dimensions. Here, labeled axes are used to give 
the dimensions of the X and Y axes and a color bar 
provides the altitude or Z axis dimensions. The soft- 
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The image of Figure 1 adjusted for altitude. 


ware operates as follows. 

The program is initiated by having the user spec- 
ify the file name of the desired image. This image is 
then transferred to the RAMTEC, and the lookup 
table is defined for a gray scale from 8 to 255. The 
user interactively indicates what portion of this image 
is to be used to create a surface. As before, software 
is used to reduce this portion to an array of *2-by-32 
points. There are a number of ways to do this. Aver- 
aging is used in this case. The portion of this image 
specified by the user is selected so that the number of 
pixels in both directions is a multiple of 32. This new 
portion of the image is then divided into 32-by-32 
sections. The pixel values in each section are totaled, 
and these totals are divided by the number of pixels 
per section. A potentia!ly undesirable effect of this 
method of reduction is low-pass filtering. High- 
frequency features in the image will be attenuated. It 
should be noted, however, that some information 
will be lost in any type of reduction. 

The program next converts the array of infrared 
data to altitude as discussed previously. The conver- 
sion table is created sometime before running this 
program and is stored in a data file which is accessed 
to make the conversion to altitude. The user specifies 
a maximum altitude that will control the size of the 
cloud surface in the vertical direction. The distance 
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from the zero elevation to the maximum altitude will 
always be the same for a given orientation. 

Since ground features are not necessarily of any 
interest and tend to create a more confusing surface, 
the program now removes all features below an alti- 
tude specified by the user by setting these array 
elements to zero. The procedure up to this point will 
produce a ragged edge on the bottom of the surface. 
To make the lower edge of the surface smooth, the 
program changes the values of the points around the 
edge of the nonzero surface to the lower altitude that 
was specified by the user. 

The program then creates the surface on the 
RAMTEC screen through the use of display sub- 
routines. The surface draws the X, Y, and Z direc- 
tions. User interaction specifies how shape-defining 
sets of contours are to be drawn. Each set of con- 
tours is chosen to be close together and cover a nar- 
row band of altitudes. The surface is drawn over and 
over again each time with a different set of contours. 
Each set of contours is drawn with a different pixel 
value which is defined as a different color. The result 
in this program is a surface on which the color 
changes from blue at low altitudes to white at high 
altitudes. This produces a known relationship be- 
tween color and altitude within the program. The en- 
tire surface, including points that are zero, is drawn 


one last time without any contours. The pixel value 
used is defined as black. The corners are made equal 


to zero to provide screen coordinates that will be used 
to draw the axes. The angle of view for the surface is 
specified by the user. 

Last of all, this program adds dimensions to the 
surface if so desired. The horizontal directions will be 
described with labeled axes. Vertical size will be in- 
dicated by a color bar with altitudes marked on it. 
The corners of the surface that were set to zero pro- 
vide the screen coordinates for the ends of the axes. 
The user must then interactively choose three of the 
four corners to define the end-points of the horizon- 
tal axes. The first corner chosen is the one with a 
screen coordinate that is lowest on the screen. The 
other two corners are those adjacent to the first. 
These are identified as right and left by comparing 
the X screen coordinates. Two lines are drawn on the 
screen for the axes using these screen coordinates. In 
order to find the scaling for each axis it must be 
determined whether the axis represents north-south 
or east-west. In this case, if the subscripts of the first 
corner are equal, then the axis on the right represents 
an east-west direction. 

The distance represented by an axis can be found 
by multiplying the distance per pixel by the number 
of pixels used to create the surface in that direction. 
The distance per pixel must be entered by the user 
and is dependent on the creation of the original im- 
age data file. Once the length represented by an axis 
is known, the interval between tick marks can be 


Figure 4. Pseudo-color representation of image in Figure 3. 
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chosen. This program chooses the interval to be some 
power of ten. This power is found by 


p = int(log,o(L)) (2) 


where 10? is the interval represented between tick 
marks and L is the total length represented by the 
axis. The distance between tick marks in terms of 
screen coordinates is found by 


AxT 


L (3) 


s = 


where s is the distance between tick marks in screen 
coordinates in the X or Y direction, A is the total 
length of the axis in screen coordinates in the X or Y 
direction, T is the distance represented by the tick 
mark interval and L is the total distance represented 
by the axis. By knowing the tick mark interval in 
screen coordinates, one can use a loop to draw tick 
marks at these intervals until a location is found for a 
tick mark that is not on the axis. The values written 
on the screen next to these tick marks are multiples of 
the interval. The value at the first corner is zero. This 
procedure is done for each axis. 

To indicate vertical dimensions, a vertical color 
bar is drawn in the upper-right corner. The color on 
this bar changes from bottom to top in the same way 
that the surface does. The relationship between the 
pixel value used to create a color on the color bar and 
its representative altitude was stored when the sur- 
face itself was drawn. The altitudes at the color 
boundaries are written on the screen. The program 
then writes the units of measurement on the screen. 

Examples of the results are shown in Figures 5 
and 6 which are three dimensional views of the two 
clouds shown in the center of Figure 4. Figure 5 is an 
intermediate result with the flat base still attached. 
Figure 6 is the finished product with some shading 
employed and dimensions included. Note that the im- 
age in Figure 6 is viewed from a different position 
than the image in Figure 5, demonstrating the option 
of different viewing angles. 

A second method for viewing clouds as three- 
dimensional objects utilizes a set of polygons to ap- 
proximate the three-dimensional surface. A software 
package called MOVIE, developed at Brigham 
Young University, is available as a starter for doing 
this. The MOVIE package contains the initial sub- 
routines needed for polygon approximation of a sur- 
face, arbitrary coloring and shading of the resulting 
surface, and the use of transparency to allow the 
viewer to see objects in back of other objects. At 
present, MOVIE does not run well on the 16 bit 
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Intermediate three-dimensional represenation of the 
two clouds in center of figure 4. 


Figure 5. 


PDP-11/60 computer because of the lack of dynamic 
memory. However, we have been able to obtain sim- 
plified approximations to demonstrate the use of the 
software for cloud representations. The dimension- 
ing software was programmed for inclusion in 
MOVIE. 

As before, the initial irnage must be sectored and 
operated on to obtain altitude information. Because 
of the lack of memory only crude polygon approx- 
imations can be generated, but these are often suffi- 
cient for cloud representations. An example of the 
results is shown for the previous two clouds in 
Figure 7. We call these slab cloud representations. 

While the polygon approximations appear to 
give a much coarser representation, the use of a 
larger computer will allow a much better polygon ap- 
proximation to be obtained. We anticipate the results 
will be more cloud-like although this hypothesis must 
be proven. 

As a final example, we have recently used our 
software to depict the output of computer models for 
modeling cloud growth. Three of the results are 
shown in Figures 8, 9 and 10 which show the begin- 
ning, an intermediate stage and the final stage of a 
simulation. 


Figure 6. Three-dimensional cloud picture with dimensions. 
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Conclusions and Future Work 


This paper describes the methods for generating 
three-dimensional cloud images from data generated 
by geosynchronous weather satellites. The resulting 
images are reasonabie representations of clouds from 
an observer’s point of view. The use of X, Y, and Z 
dimensions provides the viewer with accurate infor- 
mation concerning cloud size from the observer view- 
ing position. The cloud images can be easily manipu- 
lated to permit observation from any desired viewing 
angle and from any distance. The use of imaging to 
represent model output is also described. 

The work was done using a small computer sys- 
tem. The next phase will utilize a DEC VAX-11/780 
computer and this will greatly enhance our ability to 
investigate more sophisticated ways of utilizing and 
presenting three-dimensional data. By utilizing the 
full version of MOVIE, we will be able to more effec- 
tively represent clouds using better polygon approx- 
imations to solids and improved shading algorithms. 
The use of a transparency algorithm will allow us to 
‘*see through’’ clouds and to study cloud-cloud inter- 
actions. The transparency capability will allow the 
modeling personnel to see the inside of their model 
generated clouds, for example, stream lines can be 
displayed. 

Work will also be started to relate cloud data to 
model output. This will be done by first displaying 
real clouds and model clouds in three dimensions. 
From this we hope to learn how to use satellite data 
better for model initiaiization. We will also begin to 


Figure 7. Cloud representation using a crude polygon ap- 
proximation. 


use temporal three-dimensional cloud images to 
study methods for estimating both horizontal and 
vertical wind fields. These studies will utilize our 
present man-interactive computing facilities to allow 
Our scientists to evaluate quickly algorithmic ap- 
proaches to the problems. @ 
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Figure 9. The model cloud at intermediate time. 
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Figure 10. The model cloud at the final stage of simulation. 
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Atmospheric Boundary-Layer 
Processes Over Coastal Waters 
and Across Oceanic Fronts 


upported mainly by the Coastal Sciences Pro- 
gram of ONR and Naval Environmental Predic- 
tion Research Facility, Dr. S. A. Hsu and his col- 
leagues at the Coastal Studies Institute of Louisiana 
State University are studying the atmospheric 
changes across oceanic fronts as well as mesoscale 
processes associated with airflow adjustment over 
coastal waters. When cold air is situated over warm 
water the convective velocity scale and the height of 
the mixing layer can be estimated by the surface heat 
flux. On the other hand, when warm air flows over a 
cooler sea, the height of the modified layer is propor- 
tional to the square root of the fetch downwind from 
the coast. It can be shown that this equation is 
superior to that based on dimensional analysis. All of 
the above concepts have been verified by field experi- 
ments. At times, branches of warm oceanic currents 
such as the Gulf Stream and the Kuroshio Current, as 
well as cold oceanic eddies associated with upwelling 
such as the Somali Current, are very near the coast. 
Their impact on the atmospheric boundary layer is 
studied. It has been demonstrated that there is ap- 
proximately a 30% increase in relative humidity and 
mixing ratio from the cold side to the warm side 
about 15km downwind from the oceanic front in the 
Korea Strait. Also, a 140m increase in the mixing 
height from the cold to the warm side has been 
observed. An example is an inversion layer that ex- 
isted between approximately 200m and 600m over a 
cool pool because of the presence of a cold eddy off- 
shore of the Somalia coast in summer; no inversion 
layer was observed on the wari side of the upwelling 
zone. & 


(Dennis Conlon, ONR) 


Development of Techniques for the Remote 
Sensing of Surface Currents and Current 
Shear using Microwave Radar 


S urface currents are defined as the mean hori- 
zontal flow of water within a layer which extends 
downward from the surface to a depth of several 
meters. The accuracy of conventional moored cur- 
rent meter measurements suffers in this layer because 
of competing wave amplitude and velocity effects. 
Lagrangian drifter measurements of spar buoys pre- 
sent difficult tracking and deployment problems and 
dye studies of surface flow are difficult to interpret 
quantitatively. Still the need persists to accurately 


-map surface currents, because anything that floats on 


the surface (ships, pollutants, icebergs, nutrients and 
even waves) is transported by these currents. 
Recently, a dual-frequency microwave radar 
technique has been developed at the Naval Research 
Laboratory which is capable of remotely sensing 
both ocean surface currents and vertical profiles of 
current strength to a depth of several meters. The 
radar measures currents by monitoring the perturba- 
tions in ocean wave propagation velocity caused by 
the presence of the current. This sensitive measure- 
ment is made possible by the dual-frequency radar’s 
unique ability to ‘‘resonate’’ with one ocean wave at 
a time. The particular wave that is being used as a 
current measurement-tool can be changed at will be 
changing the separation between the pair of radar 
frequencies which are scattered from the ocean’s sur- 
face. Temporal and spatial (mapping) studies of sur- 
face currents have been carried out using this basic 
mode of operation. Modern microwave technology 
facilitates an expansion of the radar’s capability. In 
particular, signal multiplexing techniques using 
PROMS (programmable read-only-memories) have 
been developed which enable the radar to examine si- 
multaneously the effect of currents on a set of up to 
one-hundred independent ocean waves. Because long 
wavelength waves are affected by much deeper cur- 
rents than shorter wavelength waves, the multiplexed 
data can be analyzed to obtain information on the 
vertical distribution of currents beneath the ocean’s 
surface. Information about currents near the surface 
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and at depth is important because near surface cur- 
rents are ‘‘contaminated’’ by local surface winds 
which make a determination of the mean flow dif- 
ficult. 

The dual-frequency microwave current meter 
has been successfully tested from towers in the Atlan- 
tic Ocean and in the North Sea. System tests using a 
mobile platform (a research vessel the AGOR-16 
USNS Hayes) are planned for this July in both the 
Gulf Stream and in the shoals near Nantucket Island. 
(Dale L. Schuler, NRL) 


Winter-Season Thermal Stress to High- 
Latitude Reef Systems 


he effects of cold-air outbreaks over high- 

latitude reef systems are being investigated by 
scientists at the Coastal Studies Institute of Louisiana 
State University. Nan Walker, funded by a NASA 
Graduate Student Research Grant, and Drs. H. H. 
Roberts, L. J. Rouse, Jr., and O.K. Huh, funded by 
ONR, have found that these outbreaks over the reefs 
and associated shallow environments can create 
thermal stress conditions severe enough to kill reef- 
building corals, thereby affecting long-term reef 
development. Over the last five years digital thermal 
infrared data from high-resolution radiometers of 
the NOAA-5 and NOAA-6 environmental satellites 
have provided a means for determining regional sea- 
surface temperature variations during particularly 
cold outbreaks of polar air. Synoptic satellite data 
were used in conjunction with a computerized heat 
flux model to investigate shallow-water chilling 
processes. The model incorporates hourly 
meteorological parameters and the bulk aerodynamic 
equations to hindcast changes in water temperature 
throughout critical cold-air outbreaks. These results, 
coupled with in situ temperatures and biological 
observations (from Florida Bay and the northern 
Bahama Banks), have provided a data base from 
which a comprehensive appraisal of thermal stress 
events on reef communities has been formulated. 


During January i977 and 1981 consecutive cold 
fronts crossed south Florida and the northern 
Bahamas, causing abnormal chilling. Field and 
laboratory experiments identify 16°C as a thermal 
stress threshold for most reef-building corals, below 
which they have difficulty maintaining _life- 
supporting processes. Water temperatures below 
16°C were experienced in Florida Bay for 8 con- 
secutive days during the 1977 event and for 10 days in 
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1981. Minimum Florida Bay temperatures obtained 
from corrected satellite data and confirmed from 


‘‘surface truth’? measurements were 12.9°C and 
9.0°C, respectively. The figure above illustrates a 
profile of radiation temperatures extracted from 
digital satellite data along the transect A-A’. 

Massive fish kills and numerous coral mor- 
talities were observed in the study area during the 
1977 and 1981 cold-air outbreaks. Results show that 
cold-water stress conditions can exist over vast 
shallow-water areas with residence times of several 
days. Offbank movement of chilled waters, driven by 
northerly winds, density gradients, and tidal pum- 
ping, can limit coral growth and control the spatial 
distribution of reef development in high latitudes. & 


(Paul Twitchell, ONR Boston) 


Continued from page 5 


The CO, laser is focused on the water surface 
causing an explosion called ‘‘thermal blow-off.’’ 
Preliminary studies have indicated that the explosion 
creates an acoustic field in the water having a pulse 
width of 20-30 msecs. and a frequency spectrum 
peaking between 25-30 kHs. On page 5 is an example 
of a typical data return from one pulse of the radar. 
The air blast (surface explosion noise) and the 
bottom bounce signals are quite evident. The depth is 
easily calculated because the difference between the 
speed of sound in water and air is quite exaggerated. 
The measurement was taken while at sea and un- 
derway.@ 


(James S. Bailey, ONR) 61 
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